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Preface 



Quark flavor physics and lattice quantum chromo dynamics (QCD) met many years 
ago, and together have given rise to a vast number of very fruitful studies. All of 
these studies certainly cannot be reviewed within the course of these lectures. Instead 
of attempting to do so, I discuss in some detail the fascinating theoretical and phe- 
nomenological context and background behind them, and use the rich phenomenology 
of nonleptonic weak kaon decays as a template to present some key techniques and to 
show how lattice QCD can effectively help shed light on these important phenomena. 
Even though the lattice study oi K ^ tttt decays originated in the mid-eighties, it 
is still highly relevant. In particular, testing the consistency of the Standard Model 
with the beautiful experimental measurements of direct CP violation in these decays 
remains an important goal for the upcoming generation of lattice QCD practitioners. 

The course begins with an introduction to the Standard Model, viewed as an 
effective field theory. Experimental and theoretical limits on the energy scales at which 
New Physics can appear, as well as current constraints on quark flavor parameters, 
are reviewed. The role of lattice QCD in obtaining these constraints is described. A 
second section is devoted to explaining the Cabibbo-Kobayashi-Maskawa mechanism 
for quark flavor mixing and CP violation, and to detailing its most salient features. 
The third section is dedicated to the study of — >■ tttt decays. It comprises discussions 
of indirect CP violation through K^-K^ mixing, of the A/ =1/2 rule and of direct 
CP violation. It presents some of the lattice QCD tools required to describe these 
phenomena ah initio. 
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1 

Introduction and motivation 



1.1 The Standard Model as a low-energy effective field theory 

If elementary particles were massless, their fundamental interactions would be well 
described by the most general perturbatively renormalizable relativistic quantum 
field theory based on: 

• the gauge group 

SUi3)c X SU{2)l X C/(l)r , (1.1) 

where the subscript c stands for "color" , L for left-handed weak isospin and Y for 
hypercharge; 

• three families of quarks and Icptons 

u d e~ Vf. 

c s iJ~ , (1.2) 

t b T~ Vr 

with prescribed couplings to the gauge fields (i.e. in specific representations of the 
gauge groups); 

and the absence of anomalies. 

In the presence of masses for the weak gauge bosons and for the quarks and for 
the leptons, the most economical way known to keep this construction perturbatively 



renormalizable is to implement the Higgs mechanism (Englert and Brout, 1964 Higgs, 



1964), as done in the Standard Model (SM). However, this results in adding a yet 
unobserved degree of freedom to the model, the Higgs boson. 

By calling a theory renormalizable we mean that it can be used to make 'predictions 
of arbitrarily high accuracy over a very large interval of energies, ranging from zero 
to possibly infinite energy, with only a finite number of coupling constants. |^ These 
couplings are associated with operators of mass dimension less or equal to four in 3 + 1 
dimensions. 

^Beyond fixed-order perturbation theory, the U{1)y of hypercharge is trivial: the renormalized 
coupling constan t vanishes w hen the cutoff of the regularized theory is tal^en to infinity, a notion first 
suggested in I j Wilson and Kogut, 1974| . 

^If one sticks to perturbation theory, the precision reached is actually limited by the fact that 
perturbative expansions in field theory are typically asymptotic expansions. Moreover, the triviality 
of the Higgs and U{1)y sectors means that the cutoff, which we generically call A here, has to be 
kept finite. This limits the accuracy of predictions through the presence of regularization dependent 
corrections which are proportional to powers of E/A, where E is an energy typical of the process 
studied. In that sense, only asymptotically free theories can be fundamental since they are the only 
ones that can be used to describe phenomena up to arbitrarily high energies. 



2 Introduction and motivation 



Renormalizable field theories are remarkable in many ways. Consider an arbitrary 
high-energy theory described by a Lagrangian £uv (e.g. a GUT, a string theory, 
. . . ) with given low-energy spectrum and symmetries. At sufficiently low energies this 
theory is described by the unique renormalizable theory with the given spectrum and 
symmetries, whose Lagrangian we will denote £ren- Moreover, the deviations between 
the predictions of the two theories can be parametrized through a local low-energy 
effective field theory (EFT) 

rUV^Acn + EES^Of , (1.3) 

d>4 i i 

where the of ^ are operators of mass dimension d > 4 built up from fields of £ren- The 
Ai are mass scales which are much larger than the masses in the spectrum of £rcn~ 
there may be one or many of them depending on the number of distinct scales in £uv- 
The Cd^i are dimensionless coefficients whose sizes depend on how the corresponding 
operators are generated in the UV theory, e.g. at tree or loop level. 

Thus, very generally, we can write down the Lagrangian of particle physics as a 



low-energy EFT with the gauge group of Eq. (1.1), the matter content of Eq. (1.2 1 
and a Higgs mechanism: 

Cfu - Csu + j^O^Z + E E ^M'' > (1-4) 

d>6 i i 

where the left-handed neutrino Majorana mass term, ^^'^ must be 



invariant under the Standard Model gauge group (1.1 1. In Eq. (1.4), £sm is the renor- 
malizable Standard Model Lagrangian 

^SM = -Cg+f + ^flavor + ^EWSB + ■ (1-5) 

where >Cg+f contains the gauge and fermion kinetic and coupling terms, £fiavor, the 
Higgs- Yukawa terms, >Cewsb, the Higgs terms and the possible renormalizable 
neutrino mass and right-handed neutrino kinetic terms. In that sense, the renormal- 
izable Standard Model is a low-energy approximation of a more complete high-energy 
theory involving scales of New Physics much larger than M-^. 
Schematically, the gauge and fermion Lagrangian reads 

£g+f = iF;,i^r + 'A^V'. (1.6) 

It has 3 parameters, the gauge couplings (31,32153)1 a-nd is very well tested through 
experiments conducted at LEP, SLC, the Tevatron, etc. Its parameters are known to 
better than per mil accuracy. 

The Higgs- Yukawa terms are given by 

Aavor = -i^^R"^\-l/2)4'HL " + h.C. , (1.7) 

with -ij^L corresponding to the left-handed SU{2)l doublets and the right- 

handed SU{2)l singlets, associated with the — ±^ component of the doublets. 
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In this equation, (f) is the Higgs field and (f> its conjugate, (0", — 0+*). The flavor 
component of the Standard Model Lagrangian has many more couplings, 13 in fact. It 
gives rise to the 3 charged lepton masses, 6 quark masses and the quark flavor mixing 
matrix which has 3 mixing angles and 1 phase. |^ The understanding of this quark 
mixing and its associated CP violation will be the main focus of the present course. 
There is also the electroweak symmetry breaking (EWSB) contribution 

/:ewsb - - - A(0V)' • (1.8) 

It has only 2 couplings, the Higgs mass and self-coupling (/i. A), and is very poorly 
tested so far, a situation which will change radically with the LHC. 

As for the neutrino Lagrangian, little is known from experiment about its form. 
There are theoretically two possible, nonexclusive scenarios: 

1. There are no right-handed neutrinos in sight. Thus, we give our left-handed neu- 
trinos a mass without introducing a right-handed partner. In that case, Ci, = Q 



in Eq. (1.5) and we have a Majorana mass term for the left-handed neutrinos in 
Eq. (|1.4[), with 

Otl = -\^lC~^^ A^^VL + h.c. , (1.9) 



where C is the charge conjugation matrix (see Eq. (2.22 1). That is, after EWSB 
the neutrino acquire a Majorana mass through the introduction of a nonrenor- 
malizable dimension-5 operator. This implies that the Standard Model is an 
EFT and that we already have a signal for a new mass scale. Indeed, with 
nil, ~ 0.1 eV (a plausible value), eigenvalues of the coupling matrix of or- 



der 1 and ((/)) ~ 246 GeV, one finds for the mass scale M of Eq. (1.4 1 



M--^-10i^GeV (1.10) 
my 

which is tantalizingly close to a possible unification scale. 

We choose to allow right-handed neutrinos. Nr. These neutrinos must be singlets 
under the Standard Model group. Thus, they themselves may have a Majorana 
mass, but this time a renormalizable one, in addition to allowing the presence of 
a Dirac mass term: 

= NRi^Nn - {LlY^4>Nr + '^N^CM^Nr + h.c.) (1.11) 

NriI?Nr~\{vI,N^^)C ^Y^^^Jl^R^ (^1 ' ^^-^^^ 

where Ll stands for the left-handed lepton doublets, and N^^ for the charge 



conjugate oi Nr (see Eq. (2.22)). 
There are here three more possibilities: 
a) = 

In that case, the three neutrinos have Dirac masses and lepton number is 
conserved. 

^Remember that we have separated out into Ci, possible renormahzable neutrino mass terms. 
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b) M« » Y^icj,) 

Here, the see-saw mechanism comes into play: there are no right-handed neu- 
trinos in sight and aU three left-handed neutrinos acquire a mass through 
d = 5 operators: 

O - 'l'^IC4>^ (Y^iM^rX*) + O ((M«)-2) . (1.13) 
Thus, we have an explicit realization of scenario 1). 

Taking KeY,y{(j)) ~ 1 GeV in rough analogy with the r and again, ^0.1 eV, 



we obtain for the mass scale M of Eq. (1.4 1 



M ^ i^-YA^))' ^ ,Qio GeY . (1.14) 

c) Some eigenvalues of M ~ some eigenvalues of Y^, {(j)) 

Then, the sea-saw neutrino mass matrix will have more than three small 
distinct eigenvalues (actually up to six), leading to more than three light 
neutrinos. Such a possibility is constrained by phenomenology, but is not 
excluded. 

Though the topic of neutrino masses and associated mixing and CP violation is fasci- 
nating, it is not the flavor physics which is of interest to us here. Thus, this is all that 
we will say about the subject and, for the remainder of the course, we can safely take 

(5) 

— 0, forgetting about OJ^^ and C^, altogether. 

Having explored the neutrino mass Lagrangian and some of the constraints which 
neutrinos place on the scale of New Physics, we now do the same for the other com- 
ponents of /Cg^j generically denoting the scale of New Physics by A: 

1. EWSB and naturalness: besides possible right-handed neutrinos, the Higgs boson 
is the only Standard Model particle whose mass is not protected by a symmetry 
from the physics at energy scales much larger than Mw- To get a very rough 
estimate of what the contributions of New Physics to the Higgs mass could be, 
we assume that the effect of the new degrees of freedom can be approximated 
by computing Standard Model loop corrections to this mass, cutoff at a scale 
A ^ MiY that is characteristic of the new phenomena. Then, the contributions 
to the Higgs mass at one loop are given by the diagrams in Fig. |1.1| with a cutoff 
A. They yield 

SM% = j^^^iAm^,-2M^-Ml-Mfj) (1.15) 

which is dominated by the top contribution for AIh ^ 350 MeV. If, for naturalness 
reasons, we require that the physical squared Higgs mass is no less than a fraction 



/ of the correction of Eq. (1.15), then we find that 



700 GeV (Afe/115GcV) ,^ 

(1.16) 



V7 



A/ff-115 GeV ^ 
310 GeV ) 
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70 M/± 



I 



t 

Fig. 1.1 Diagrams wliicli contribute radiative corrections, 5M%, to the Higgs mass squared 
at one loop 



Thus, if we allow at most 1% of fine tuning on the Higgs mass squared, Eq. (1.16 1 
says that new physics must appear below A„at ~ 7TeV. 

Gauge sector and flavor conserving d = 6 operators: consider, for instance, Owb = 
gig2{<l>^ cr°'(f)W^^B^i,, which couples the W bosons to the U{1)y gauge boson B. 
Precision electroweak data, assuming that the Wilson coefficient of this operator is 



of order one, impose the following constraint on the scale of New Physics ( Barbieri 



et al., 2004) 



A>5TeV 95% CL. (1.17) 

3. Flavor physics and, in particular, flavor changing neutral currents (FCNC): con- 
sider, for instance, K^-K^ mixing. In the absence of electroweak interactions, the 
long-lived is a CP odd state, whereas the short-lived is a CP even state. 
When these interactions are turned on, these two degenerate particles acquire a 
minuscule mass difference, which is measured experimentally to be: 

AAfx = Mk,^ - Mks ^ 3.5 X 10-^^ MeV . (1.18) 

Consider now the contribution to AMk of an arbitrary d = 6, AS = 2 operator 
schematically written as {ds){ds): 

AMk = 2 X ^^'=^ff^{K\-sd){mK') . (1.19) 

Now, assuming that R.eC^5^2 ^ 1 that the matrix element is of order the 
fourth power of a typical QCD scale, e.g. ^ M^, we get 

A/2 

A > ^ ~ 10^ TeV (1.20) 

^/Mk AMk 

which is orders of magnitudes larger than the lower bound imposed by gauge 
sector and flavor conserving transitions, as well as than the upper bound imposed 
by naturalness. 

Thus, if we do not make any assumptions about how the New Physics breaks flavor 
symmetries, we are forced to push this physics to very high scales. Said differently, 
flavor physics is sensitive to very high energy scales if the New Physics is allowed to 
have a flavor structure which differs from that of the Standard Model. Therefore, one 
assumption commonly made is that the New Physics breaks the flavor symmetries 
with the same Yukawa couplings as in the Standard Model. This assumption is called 
Minimal Flavor Violation (MFV). For instance, we might have, in the case of K^-K'^ 
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mixing, the following operators contributing: ^(sflXsddi)^, j2{sLY*jYsdJp.dL)'^ , etc. 
I will leave you work out the corresponding scales, A, but they are certainly much 
lower and in line with those obtained from flavor-conserving physics. 



1.2 Flavor physics phenomenology 

As we shall see shortly in more detail, the Standard Model has a very rich and con- 
strained flavor structure, which includes: 

• mixing of quark flavors; 



• CP violation by a unique invariant J (Jarlskog, 1985), discussed in Sec. 2.3.4 

• the absence of tree-level flavor changing neutral currents (FCNC) . 
All of these features are encapsulated in: 



the Cabibbo-Kobayashi-Maskawa (CKM) matrix (Cabibbo, 1963 Kobayashi and 



Maskawa, 1973) 



(Vud Vus Vub\ 
Vcd Vcs Vcb , (1.21) 
Vtd Vts Vtb ) 

which is unitary. It has 3 mixing angles and a single phase, which is responsible 
for CP violation. 
• the quark masses: vriq, with q — u,d, s, c, 6, t. 

For their discovery, in 1973, that Nature's CP violation and rich flavor structure can 
be well described when a third generation is added to the SU {2)lxU{1)y electroweak 



model (Kobayashi and Maskawa, 1973), Kobayashi and Maskawa were awarded part 
of the 2008 Physics Nobel Prize. 

Because this flavor structure is so intriguing and most probably contains important 
information about physics at much higher energies than currently explored, particle 
physicists have invested a considerable amount of effort in exploring it theoretically 
and experimentally over the past five decades. This exploration has multiple goals: 

1) To determine from experiment the matrix elements of the CKM matrix V , which 
are important parameters of our fundamental theory. 

2) To verify that the CKM description of quark flavor mixing and CP violation is 
correct, e.g.: 

• Can all of the observed CP violation in the quark sector be explained in terms of 
a single phase? 

• Is the measured matrix V unitary? 

The latter can be tested by verifying whether 

\yuD?^l and \VuD? = l. (1.22) 

D—d.sM U—u,c± 



If either of these sums turn out to be less than 1, that would signal an additional 
generation or family. On the other hand, if either one is larger than 1, completely 
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Fig. 1.2 The (db) unitarity triangle. 



Fig. 1.3 From top to bottom, the (db), (sb) and (ds) unitary triangles, normalized by 
V'cDi '14"d2/I^'=Oi K-D^ 1 with (Di,Z)2) = {d,b), {s,b) and (d, s) respectively, and drawn to a 
common scale. 

new physics would have to be invoked. The unitarity of V also implies that the scalar 
product of any two distinct columns or rows of the matrix must vanish, i.e. 

E ^UD.V^D, = 0, for D,^D2, (1.23) 

U=u,c,t 

E VuroV^.D = 0, for Ui^U2. (1.24) 

D=d,s,b 

These relations can be represented as triangles in the complex plane, which are tradi- 



tionally labeled by their unsummed flavors, i.e. (Di,D2) for those of Eq. (1.23) and 



{Ui,U2) for those of Eq. (1.24 1. In the absence of CP violation, unitarity triangles 
would become degenerate. The (db) triangle is shown in Fig. 1.2 In Fig. 1.3 it is 
drawn to scale with two other triangles to give you a sense of the variety of unitarity 
triangles and the difficulties there may be in measuring some of their sides and angles. 

The (db) triangle has been the focus of considerable experimental (LEP, S-factories, 
Tevatron, . . .) and theoretical (QCD factorization, Soft Collinear Effective Theory 
(SCET), lattice QCD (LQCD), . . .) effort in the last ten to fifteen years. With the 
arrival of the LHC and, in particular, the experiment LHCb, the focus is shifting from 
the study of the towards the study of the Bg meson and thus toward investiga- 
tions of the (sb) triangle. Here too, lattice QCD has a considerable role to play, most 
notably in the study of the Bg-Bg mass and width differences, of the leptonic decay 
Bg — !■ IJ.~^ fJ-~ or of the semileptonic decay Bg — )■ (p^^ . 

The strategy here is to verify the unitarity of the CKM matrix by performing 
redundant measurements of: 
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• triangle sides with CP conserving decays, 

• angles with CP violating processes, 
and checking that the triangles indeed close. 

3) To determine in what processes there is still room for significant New Physics 
contributions. For instance, 0(40%) effects are still possible in B^-B^ mixing from 



New Physics with a generic weak phase (Lenz et ai, 20101. 



4) To constrain the flavor sectors of beyond the Standard Model (BSM) candidates. 
As we saw above, it is difficult to add new physics to the Standard Model without 
running into serious problems in the flavor sector. 

5) To actually find evidence for beyond the Standard Model physics. This is most 
likely to be found in processes which are highly suppressed in the Standard Model, 
such as FCNC. 

6) If new particles and interactions are discovered, it is important to investigate their 
quark and flavor structure. 

All of these goals require being able to compute reliably and precisely flavor ob- 
servables in the Standard Model or beyond. A high level of precision has been reached 



already on the magnitudes of individual CKM matrix elements (Charles et al., 2005) 



/0.97425t!^ 
0.22529 



0.00018 



t 



0.00077 
0.00077 



0.22543 
0.97342 



V0.00858t|] °E;°i^° 0.04054+^ r;?a 0.999141 



s 

+0.00077 
0.00077 

+0.00021 
-0.00019 

^+0.00057 
*-0. 00129 



0.00354: 
0.04128 



+0.00016 
-0.00014 



+0.00058 
'0.00129 



+0.000053 
-0.000024 



(1.25) 



assuming the correctness of the Standard Model and, in particular, CKM unitarity. The 
most poorly known CKM matrix elements are \ Vub\ and \Vtd\, both with an uncertainty 
around 4%. Then come |Vts| and \Vcb\, with an uncertainty of about 2%. Thus, to have 
an impact in testing the CKM paradigm of quark flavor mixing and CP violation, and 
to take full advantage of LHCb results, the precision of theoretical predictions must 
be of order a few percent (better in many cases). This is no small challenge when 
nonperturbative QCD dynamics is involved. 



1.3 Flavor physics and lattice QCD 

Lattice QCD plays and will continue to play a very important role in flavor physics, 
by providing reliable calculations of nonperturbative strong interaction corrections to 
weak processes involving quarks. 

The processes for which LQCD gives the most reliable predictions are those which 
involve a single hadron that is stable against strong interaction decay in the initial 
state and, at most, one stable hadron in the final state. Resonances (i.e. unstable 
hadrons) are much more difficult to contend with, especially if many decay channels 
are open. Similarly, final states with more than a single stable hadron are much more 
difficult, especially if these hadrons can rescatter inelastically. This will be discussed 
in much more detail in Sec. 13.71 
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Now, to determine the unique CKM matrix phase or, more precisely, the CP vio- 
lating parameter J, lattice QCD can have an important impact through the following 
processes: 
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• Indirect CP violation in if — tttt decays. This occurs through the process of 
K'^-K^ mixing, which is given by the imaginary part 

Im 



where the eUipsis stands for the other box contribution. It is a \/S.S\ = 2 FCNC 
that will be discussed in detail below. 
• Direct CP violation in K ^ tttt decays. The processes which contribute are given 
by the following \/S.S\ — 1 amplitudes 

ir+ rf 



K" I 11, c. t 

rf 



+ 




in the case of K'^ — > tt+tt^, where the ellipsis stands for missing diagrams similar 
to those drawn. The diagrams for K'^ — tt^tt" are analogous. Again, this particular 
phenomenon will be discussed in detail in the sequel. 






1.4 Low-energy effective field theories of the Standard Model 

With present knowledge and present computer resources, we cannot simulate the full 
Standard Model in lattice field theory calculations. In particular various degrees of 
freedom must be "eliminated" from the calculations for the following reasons: 

• W, Z and t: there is no hope to be able to simulate these degrees of freedom whose 
masses are Mw,z ^ 80 90GeV and nit ~ 175 GeV on lattices which must be 
large enough to accommodate 135 MeV pions, i.e. with sizes L > A/IVU - 6fm.[^ 
Since we would also have to have arrit <C 1, with a the lattice spacing, to guarantee 

*The factor of 4 in is a conservative rule-of-thumb estimate which guarantees that finite- 

volume corrections to stable hadron masses, proportional to e~^^'^^, are typically below the percent 
level. 
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Fig. 1.4 Examples of correlation functions required for the lattice computation of the K^-K'^ 
mixing amplitude. The diagram on the left exhibits the type of three-point function with a 
four-quark operator insertion, required to obtain the K'^-K'^ amplitude when the charm quark 
is integrated out. The diagram on the right shows a four-point function which is required when 
the charm quark is kept active. 

controlled discretization errors, the number of points on the lattice would have 
to be L/a ^ 4mt/M^ ~ 5.2 x lO'^, which is beyond any foreseeable computing 
capabilities. Perhaps even more important, however, is the fact that we just do 
not know how to discretize nonabelian, chiral gauge theories (please see David's 
lecture notes in this volume (Kaplan, 2010)). 



h: even the b quark, with rub ^ 4.2 GeV, would require lattices with L/a 120, 
which is already too much for present technology. 

c: with TOc ^ 1.3 GcV, the charm is a borderline case, both on the lattice and 
in terms of QCD. On the lattice it can be included in simulations, but with 
arric ~ 0.35 at best, discretization errors remain an important preoccupation. 
From the point of view of QCD, the charm is not quite a heavy quark-heavy 
quark effective theory is only marginally applicable since rric is not much larger 
than typical QCD scales-and it is clearly not light-it is not in the regime of chiral 
perturbation theory. So its inclusion should be considered with care. In addition, 
its inclusion in weak processes can mean significantly more complicated correlation 
functions to compute. For instance. Fig. 1.4 illustrates the type of correlation 
functions required to determine the amplitude for K'^-K^ mixing, in the absence 
of charm (diagram on left) and in the presence of a dynamical charm quark 
(diagram on right). In the absence of charm, it is a rather standard three-point 
function which must be computed while in its presence, it is a four-point function, 
with two four-quark operator insertions. However, in certain circumstances, the 
inclusion of a dynamical charm quark significantly simplifies the renormalization 
of the weak effective theory, as briefly discussed in Sec. 3.8 for the A/ =1/2 rule. 



Fortunately, with the precisions required at present and in any foreseeable future, 
it is not necessary to include virtual W , Z, t and h contributions. The situation with 
the charm is less clear, as {Mp/rricY sea corrections can, in principle, play a role when 
percent precisions are reached. So in considering processes involving these massive 
particles, we can turn to effective field theories (EFTs) in which the W ^ t, b and 
possibly c and r are no longer dynamical degrees of freedom. Thus, we are left with an 
SU{3)c X [/(1)em gauge theory of color and electromagnetism. This theory includes 
the virtual effects of the following degrees of freedom: 
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• 2 massless gauge bosons: the qluon, g, and the photon, 7 

• 3 to 4 quarks: the u, d, s and possibly the c 

• 2 to 3 leptons: the e, /x and possibly the r 

• 3 neutrinos: the Vg, v^^ and the Vt- 

It also includes local operators of naive mass-dimension d > 5, which result from 
integrating out the heavy W, Z and t. Moreover, in this theory, the 6-and possibly the 
c and the r-arc described by heavy fcrmion effective theories, in which the antiparticle 
of the fermion considered is integrated out. It is at this level that LQCD enters, to 
describe the nonperturbative effects of the strong interaction. 



2 

Standard Model and quark flavor 
mixing 



As already discussed, the Standard Model has a highly constrained quark flavor struc- 
ture, parametrized by the CKM matrix and the quark masses. We will now see more 
precisely how it arises and what its basic implications are. 

2.1 On the origin of quark flavor mixing in the Standard iVIodei 

Let us look in more detail at the quark and gauge sectors of the Standard Model. With 

the notation (dim. rep. SU{3)c, dim. rep. SU{2)]^)y, where "rep." stands for represen- 
tation, the quark content of the Standard Model is given by 

Ur ~ (3, 1)2/3 (2.1) 
£>fl~ (3,l)_i/3 , 

with U = {u,c,t)'^ and D = {d,s,b)^. The coupling of these quarks with the gauge 
bosons is given by 

£g+q = • • • + QlPQl + UrPUr + DrPDr , (2.2) 

where the ellipsis stands for the gauge kinetic and self-interaction terms. This La- 
grangian has the following global symmetries: U{3)l on Ql,U (3)f/„ on Ur and U (3) Dr 
on Dr. 

Now let us look at the quark Yukawa terms. After spontaneous symmetry breaking 
by the Higgs field, keeping only the terms proportional to the Higgs v.e.v. (0), we 
have: 

CI -UrMuUl - DrMdDl + h.c. , (2.3) 

where Mu and M75 are arbitrary 3x3 complex matrices. On these terms we perform 
the following set of flavor transformations which leave jCg+q invariant: 

• the SU{3)u„ transformation: Ur — V^Ur, 

• the SU{3)l transformation: Ul ^ V^Ul and L»l V^D^, 

• and the SU{3)]jf^ transformation: Dr V^Dr, 
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with y^, Vi; and such that: 

= V^^MuVE , (2.4) 

and 

V^^Md = M^V^^ (2.5) 

where Mjj are diagonal matrices with real positive entries which are the quark 
masses. The scx-ond (equation defines a fourth unitary matrix . Under these ro- 
tations, the quark mass Lagrangian transforms as 

£1 ^ -UrM^Ul - DrMI,[V^WEDl] + h.c. . (2.6) 

The up quark mass matrix is diagonal in this basis, but not the down quark matrix. 
In addition, we have exhausted the flavor transformations allowed by £g+q. Thus, if 
we want to work in a mass basis (i.e. a basis in which all quark masses are diagonal), 
we have to perform the additional transformation: 

• Dl ^ vPv[>Dl, 

which is not a symmetry of £.g+q. Clearly, the only terms in £g+q which are affected 
are those which couple Ul and Dj^ . They are transformed in the following way: 

^UL\f^+^VDL + h.c. , (2.7) 

v2 

where CC stands for charged current and V = V^^V^ is the CKM matrix. All other 
terms are left unchanged. In particular, the neutral currents: 

+ |(t7i,5z.)^(^^) (2.8) 

+ '^giUR^UR - ^q^Dr^Dr , 

which are flavor diagonal, remain diagonal in the mass basis. Moreover, the Higgs 
couples to fermions through their masses and therefore has a diagonal coupling to the 
quarks in the mass basis. The fact that all uncharged couplings remain diagonal in 
the mass basis implies that there are no tree level FCNC transitions in the Standard 
Model. 

2.2 Properties of the CKM matrix 

In this section, we look in more detail at what are the key properties of the CKM 
matrix. 



Properties of the CKM matrix 15 















b' 







2.2.1 Degrees of freedom of the CKM matrix 

The flavor eigenstates, {d' , s' , b'), are related to the mass eigenstates, {d, s, b), through: 



(2.9) 



where the CKM matrix is unitary, i.e. V^V = VV^ = 1. Since F is a 3 x 3 complex ma- 
trix, it has 9 phases and 9 moduli. Unitarity imposes 3 real and 3 complex constraints. 
Thus, we are left with 6 phases and 3 moduli which, because of the normalization of 
the rows and columns to one, can be written as angles. 

Now, aside from the CC interactions, all other terms involving quarks are diagonal 
in flavor, and have LL, LR, RL and RR chiral structures. Thus, we can perform vector 
(not axial) phase rotations on each flavor and leave all of these other terms invariant. 
However, under the phase rotations 

UL,R^e'^"UL,R and Dl,r ^ e'^" Dl,r , (2.10) 

with U = u, cov t and D = d, s ov b, the U D component of the CKM matrix transforms 
as 

VuD^VuDe'^^''-^''K (2.11) 

These transformations can be used to eliminate phases in V. Although there are 6 
phases 9jj and Od, only 5 phase differences 6d — Ojj are independent. Thus, only 5 of 
the 6 phase can be eliminated. 

We have now exhausted the field transformations that can be used to reduce the 
CKM matrix' degrees of freedom. Thus, V has 3 angles and 1 phase. 

2.2.2 Standard parametrization of the CKM matrix 

The idea behind this parametrization is to write V as a product of 3 rotations between 
pairs of generations, throwing the phase into the 1^3 rotation, so that it multiplies 
the smallest mixing coefficients. Thus, 

V = R32 diag{l, 1, e'^}i?3i{l, 1, e-"}i?2i , (2.12) 

with the rotations 

(C12 S12 o\ 
-S12 C12 , (2.13) 
ij 

C12 = cos ^12 and S12 = sin ^12, and similarly for the other rotations. This yields the 
following expression for the CKM matrix: 

C12C13 C13S12 sise" 

V = I -S12C23 - Ci2S23Sl3e"' C13C23 - Sl2S23Sl3e*'' C13S23 | , (2-14) 
S12S23 - Ci2C23Ci3e*'' -C12S23 - C23Sl3Sl3e*'' C13C23 



and the angles are chosen to lie in the first quadrant. Note that this parametrization 
is not rephasing invariant. 
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2.2.3 Wolfenstein parametrization 

Experimentally, it is found that 1 >• S12 >• S23 ^ Si3, i-S. mixing gets smaller as one 
moves off the diagonal. It is convenient to exhibit this hierarchy by expanding V in 



powers of S12, i.e. in the sine of the Cabibbo angle 612 (Cabibbo, 1963 Gell-Mann 



and Levy, 1960). This yields the Wolfenstein parametrization (Wolfenstein, 19831. To 



implement this expansion, we define (Buras et al., 19941 

A = S12 



Wu. 



AX-" 



S23 



cb\ 



(2.15) 



uh 



and make the appropriate replacement in the standard parametrization of Eq. ( 2.14 1. Then, 



V 



1-AV2 A AX^ip-ij-iY 

-X 1 - AV2 AX^ 

AX^{1~ p~iri) ~AX^ 1 



(2.16) 



which clearly exhibits the hierarchy of mixing. 
2.2.4 CP violation 

Let us now see how CP violation arises in the Standard Model. Under parity, the 
charged W^-boson fields transform as 



with xp = (.t",— a;). Similarly, under charge conjugation, 
Thus, under CP, these bosons transform as 



(2.17) 



(2.18) 



(2.19) 



In my favorite Dirac spinor basis, the parity transform of a fcrmion field is given by: 



In this basis, 



r 



■>pR 



o-'' 
a-^ 



(x) 



{xp) 



and 



7^ = 



-I 
/ 



(2.20) 



(2.21) 



where / is the two-by-two unit matrix, cr'^ = (/, f?), cr^ = {I,— a) and a are the 
Pauli matrices. Clearly, parity is not a symmetry of the Standard Model since left and 
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right fermions belong to different representations of the Standard Model group. Under 
charge conjugation, we have 

^{x) = J;^) (x) A = [lfj^*J (^) = ^7^V^(^) = C^P^ix) . (2.22) 

Again, charge conjugation is clearly not a Standard Model symmetry. However, the 
CP operation, 

has a chance of being a symmetry transformation as it does not mix left and right- 
handed fields. Using the well known Pauli matrix identity, cr^crV^ = —a" and the 
anticommutation of fermion fields, the CC quark term of the Standard Model La- 
grangian transforms, under CP, as: 

(2.24) 

Since V* ^ V in the presence of a nonvanishing phase, S, CP is potentially violated 
in the Standard Model. We will see below what the necessary conditions for CP to be 
violated are. 



2.3 CP violation and rephasing invariants 



The standard parametrization of the CKM matrix V, given in Eq. (2.14), corresponds 
to a particular choice of quark field phases. Observables cannot depend on such choices. 
Therefore, it is important to find rephasing invariant combinations of CKM matrix 
elements. 

2.3.1 Quadratic invariants 

The moduli 

luh ^ \Vud\' , (2.25) 

with U = u,c,t and D ~ d,s, b, are clearly rephasing invariant. There are 9 of these. 
Now, unitarity requires that: 



= Y.U=u,c,t ^UD = 1 



^='^^'''> ; , (2.26) 



(2) 

which yields 6 constraints on the However, we clearly have J2u '^u = ^d'^d, 

which means that there are only 5 independent constraints. In turn, this means that 

(2) 

there are 4 independent quadratic invariants Ijjj^, which are obviously real. 
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2.3.2 Quartic invariants 

We now define 

^u!diU2D2 = ^UiDiVu2D2VuiD2^U2Di ■ (2.27) 

These products of CKM matrix elements are also clearly rephasing invariant, since for 
every field which one of the CKM factors in (2.271 multiplies, another factor multiplies 
its Dirac conjugate. In Eq. (2.27), Ui, U2 {Di, D2) are chosen cyclically amongst u, c, t 
{d,s,b) so as to avoid I^^-* = (I^^^)^ as well as to avoid obtaining complex conjugate 

invariants, e.g. ItjId2U2Di ~ -^tjId^U2D2- With these constraints, there are 9 invariants. 
However, not all of these invariants arc independent. Indeed, unitarity yields 

T.D=d.s,b'^UiDV{j^jj = Ui^U2 , , 

j:u=u,c,tVuD,V{}^,=0 D,^D2 ■ ^ • 

This implies, in turn: 

yUiDiV{}^;[j^ ^Y.U^Ui^UDiV{}jj^ (II) 

Multiplying both sides of (I) by Vu^D^Vfj^j^^ while maintaining the cyclicity of indices, 
yields 

^ulot U2D2 ^ ~\^UiDt Vu2 D2? ~ Iu!d2 U2 D3 ■ 
Similarly, multiplying both sides of (II) by VU2D2VU2D1 gives 

Iu!d,U2D2 = -\VU2D,VU2D2 P - I^%,U3D2 (IV) . 

Thus, the 9 /'"^^'s can all be written in terms of I^^J^g, for instance, and the 4 inde- 
pendent ijjl)- Moreover, (III) and (IV) imply that all 9 I'^*-' have the same imaginary 
part. 

2.3.3 Higher-order invariants 

Higher-order invariants can, in fact, be written in terms of J^^^s and J^^^s. For in- 
stance, the sextic rephasing invariant Vc/jnjT^;72D2^c/3D3^(7i£)2^c/2-D3^[/3Di equal to 
^UiDiU2D2^U2DiU3D3/ ^UiDi - "^^^s obviously fails in singular cases, but these will not 
be considered here because they are irrelevant in practice. 

2.3.4 Jarlskog's invariant 

This whole discussion of invariants implies that there is a unique, imaginary rephasing 
invariant combination of CKM matrix elements. This invariant must appear in all 
CP violating observables, because it is the imaginary component of the CKM matrix 
which is responsible for this violation. This invariant is known as the Jarlskog invariant 
(Jarlskog, 1985|: 



J = Im/*") = cl^Ci2C23Si2S23Si3Ss 

^ X'^A^r^ + OiX^^) , (2.29) 

where fj = 77(1 — A^/2). This, in turn, means that to have CP violation in the Standard 
Model, 6112, O23, must not be or 7r/2 and d ^ 0,Tr. Moreover, CP violation is 
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maximal for ss = 1, 0i2 = ^23 = 7r/4 and S13 = 1/a/3- At that point, the Jarlskog 
invariant takes the value 

Jrna. = ^ 0.1 . (2.30) 



However, global CKM fits yield (Charles et al, 2005) 



J = 2.96^1? X 10"^ < J„ 



(2.31) 



i.e. CP violation in Nature is very far from being as large as it could be. 

Instead of looking at the CKM matrix for a rephasing invariant measure of CP vio- 



lation, we can go back to the mass matrices Mu and Md- In Sec. 2.1 after performing 
only flavor transformations which are symmetries of Cg^q, we reached the point where 
the quark mass term was given by Eq. (2.6) 



CI, ^ -UrM^Ul - DrM^V^Dl + h.c. , (2.32) 
where V is the CKM matrix. Now, by performing the flavor symmetry transformation 

Dr V^Dr , (2.33) 



we obtain 



CI, ^ -UrM^Ul - DrVM^V^Dl + h.c. 



(2.34) 



In this equation, both M^j and Afj^ = VMfjV^ are hermitian matrices. This means 
that the commutator Cj = [M^, M'j^] is pure imaginary. Since the only source of imag- 
inary numbers in Eq. ( 2.34[ ) is the phase of the CKM matrix, Cj carries information 
about CP violation in the Standard Model. However, Cj is not a rephasing invariant. 



Defining the matrices Pjj = diag{e 



under the rephasing operations of Eq. (2.10), Cj transforms as 



'} and Pd = diag{f 



Cj 



(2.35) 



where Pjj cancels against P^j in the first term, because and Pu are diagonal. 
Nonetheless, det Cj is rephasing invariant, because and Po are diagonal. Thus, 
following Jarlskog ( Jarlskog, 1985[ ), we consider 



dot Cj ~ 2iJ X Y[ 

Ul>U2 



I n ("^^1 

Di>D2 



(2.36) 



In light of what was discussed in the section on rephasing invariants, det Cj must 
be the only imaginary, rephasing-invariant quantity that can be obtained from the 



mass term of Eq. (2.34), which is the most general mass term that can be written for 



quarks in the Standard Model. This means that the presence of CP violation is the 
Standard Model is equivalent to det Cj ^ 0. In turn this implies that there will be 
CP violation if and only if the conditions on the mixing angles and the phase of the 



CKM matrix given after Eq. (2.29) are obeyed, but also if and only if there are no 
mass degeneracies in the up and down quark sectors. 
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2.3.5 Unitarity triangle areas 

As we saw earlier, the unitarity of the CKM matrix gives rise to 6 unitarity triangles, 
defined by Eqs. ( 1.23 )-( 1.24). The areas of the {Di,D2) triangles are given by 



l\-^lmiVuD,V:D.V:o,VcD.)\ 
-Im/^) - - J 



(2.37) 



where the last line follows from Eq. ( 2.29 1. Similarly the areas of the {Ui, U2) triangles 
are 



Thus, all 6 triangles have the same area, which is given by the Jarlskog invariant. 
Since CP violation can only arise if J 0, none of the unitarity triangles can be 
degenerate if Standard Model CP violation is measured in Nature. 



3 

A lattice case study: K tvtv, CP 
violation and A/ = 1/2 rule 



Having introduced the Standard Model and flavor physics, we now turn to an impor- 
tant set of processes that have been nagging theorists for over four decades: K nn 
decays. These decays have been a rich source of information and of constraints on the 
weak interaction. In 1964 they provided the first evidence in Nature for indirect CP 
violation, which arises in the mixing of the neutral kaon with its antiparticle before 



the decay into two pion (Christenson et ai, 1964). Then, in 1999, CP violation which 



arises directly in the flavor-changing decay vertex was discovered in these same decays, 



after more than 20 years of experimental effort (Fanti et a/., 1999 Alavi-Harati et al, 



1999) 



Currently these decays still give very important constraints on the CKM paradigm, 
through the measurement of indirect CP violation, parametrized by e. And as far 
as we presently know, direct CP violation in these decays, parametrized by e', may 
be harboring New Physics. Moreover, K — >■ tttt decays display what we believe are 
unusually large, and certainly poorly understood, nonperturbative QCD corrections, 
which go under the name of AI =1/2 rule. 

At first sight the study of these decays is a perfect problem for the lattice. Only 
u, d and s (valence) quarks are involved, so that one expects controllable discretiza- 
tion errors. Of course, pions are light, which makes them diflicult to simulate, and 
there are also two hadrons in the flnal state. But SU{i) chiral perturbation theory 
(xPT) at LO relates K ^ tttt to K tt and K {) amplitudes which are simpler to 



compute (Bernard et a/., 1985). Moreover, xPT at NLO relates X — tttt ampfitudes 



obtained with heavier pions to the same amplitudes with physically light pions ( Kam- 



bor et al., 1990 Kambor et al., 1991 1. Given the typical size of chiral corrections, one 
would expect that we could at least be able to get an 0(20 — 30%) estimate of the 
relevant amplitudes. In addition to which we might expect good signals since only 
pseudoscalar mesons are involved. 

Despite all of the positive indications that the lattice should be able to provide 
valuable information about these decays, all attempts to account for nonperturbative 
strong interaction effects have failed, except in the study of indirect CP violation. 
And though much progress has been made on many aspects of these decays over the 
years, providing a fully quantitative description still remains an open problem. Thus, 

'^This was certainl y true at the time of the school, but the sit uation has been moving quite fast 
since then [please see i jChrist, 2010a||Liu, 20To| [Sachrajda, 2010a| for an update]. 
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I have chosen to focus on this particular topic in my discussion of the appHcation of 
lattice methods to flavor phenomenology. 

3.1 K ^ nn phenomenology 

Kaons have strong isospin 1/2, while pions have isospin 1. The weak decays of a kaon 
into two pions can thus occur through two channels in the isospin limit: 

• the A/ = 3/2 channel, where the final two pions are in a state of isospin / = 2, a 
state which we label (tttt) 2; 

• the A/ = 1/2 channel, where the final two pions are in a state of isospin / = 0, a 
state which we label (7r7r)Q. 

Decay into a two-pion state with isospin / = 1 is forbidden by Bose symmetry. 
We denote the amplitudes for K tttt decays by: 

T[X" ^ (tttt)/] = lAie'^' , (3.1) 

where Si is the strong scattering phase of two pions in the isospin /, angular momentum 
J — channel, defined through 

T[(7r7r) / ^ (tttt) /] = 2e'^' sin 61 . (3.2) 



In Eq. (3.11, is the flavor eigenstate with I3 = —1/2 and strangeness S' = 1: it is 
composed of a d and an s quark. 

Using this notation, we have the following isospin decompositions for the K ^ tttt 
amplitudes: 

- iT[K" ^ TT+TT-] = ^A^e^^ + -^^oe'*" 
\/6 \/3 

-iT{K'' ^ ^V"] = \/ ^^26'*^ - ^^oe"^" (3.3) 
V 3 \/3 

-iT{K^ ^ tt+tt"] = ^A^e^^ , 

where the coefficients of the various amplitudes are simply SU{2) Clebsch-Gordan 
coefficients. If CP violation is present, then A*j Aj. 

Now, in the absence of CP violation, the two physical neutral kaon states Ks and 
Kl are also CP eigenstates: |^ 

\Ks/l) ^ \K^) ^ ^ T \K')) , (3.4) 

with CP\K±) = ±\K±). The CP even Ks decays only into two pions, while the CP 
odd decays into three. Because of the phase space available to the decay products, 

^Here and in the following, we will assume that CPT is conserved. We will also work in the strong 
isospin symmetry limit. 

^The CP transformation of the neutral kaon states is chosen here to be CP\K^(p)) = — \K^{pp)) 
and CP\K^{p)) = —\K'\pp)) , where pp is the parity transformed four-momentum p = {p', —p). In 
the following, we will ignore the momentum labels unless they play a relevant role. 
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the former is much shorter hved than the latter, with a hfetime rg ~ 10~^° s versus 
tl ^ 5 • 10^^ s. This explains the subscript S for short and L for long. 

In Nature, the weak interaction breaks CP, and Ks and Kl are not pure CP 
eigenstates. As a result of K'^-K'^ mixing through the weak interaction, we have 

I^V5> = J-j^ (|i^T> + ~^\K±)) (3.5) 
((l + ?)|X")±(l-?)|i?0)) , (3.6) 



V2./T+W 

where e is a small complex parameter. 

The neutral kaons form a two state quantum mechanical system which can be 
described by a nonhermitian, 2-by-2 Hamiltonian: 

Hij = M^j — -Tij , (3.7) 

where i,j = l corresponds to K'^ and i,j — 2 to K'^. CPT implies that Hu = H22 and 
H21 = Hi2- To determine the elements of this matrix, we first decompose the effective 
Hamiltonian for the Standard Model into a QCD+QED part, i?QCD+QED7 and a weak 
part, and work to second order in the weak interaction. Then, 

where M^q is the mass common to and A"", as given in QCD and QED, and En 
is the energy of the intermediate state \n). 

Now, Cauchy's theorem implies (with P the principal part) 

P ( ^—-] ~ inS{E - oj) , (3.9) 



uj — E + ie — E ^ 

where the first term on the RHS of this equation will yield the dispersive contribution 
to the Hamiltonian of Eq. ( 3.7 1 (i.e. the mass term) and the second term, the absorptive 
part (i.e. the width term). Then, the off-diagonal element of the mass matrix is 

{K'^\Has=2\K°) , 1 „v/^ {K"\H^s=i\n){n\H^s=,\K'') 



M12 = 



2Mko 2M 



Iko MKO-En ' ^ ' ' 



where the term with the double insertion of the AS — 1 Hamiltonian gives rise to 
long distance contributions, since the states \n) which can contribute are light. For 
instance, \n) can be a tt+tt^ state. 

The off-diagonal element of the width matrix is given by the absorptive part of the 
integral: 

ri2 = 7^i {K"\HAs=i\n){n\HAs=i\K''){2n)S{En - M^o) (3.11) 
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Now, the physical states Kj^ and Ks are the eigenstates of Hij with eigenvalues 
Ml — fFi and Ms — fTs, respectively. It is straightforward to express these quantities 



in terms of the and Tij . 



AM, 



K 



M 



K, 



one obtains: 



Defining 

- Mks 

1 + e , 



and 



Mr 



Ar.- = r 



Ks 



1 - e 



AAf, 



if 



;Ar 



(3.12) 



(3.13) 



K 



Solving for e gives us explicit expressions for computing the relationship of the physical 
eigenstates Kl and Ks to the flavor eigenstates and through Eq. (3.6). It is 
worth noting that 



AM 



K 



2Mi 



and 



Ar 



K 



(3.14) 



to first nontrivial order in the weak phases. 

Having implicitly worked out the relation between physical and flavor eigenstates, 
we return to i^T — > tttt decays. What is actually measured are the amplitude ratios 



?7oo 



_ T[Kl 
T\Ks 



and 



_ T\Kl 
T[Ks 



(3.15) 



Experimentally, |?7ool — 2 x 10 ^, and |?7oo/^h I — 1 (Nakamura et al., 2010). These 



ratios are clearly CP violating since i^L does not decay into two pions if CP is con- 



served. The CP violating decays Kl tttt can occur in two ways. As seen in Eq. (3.5), 
\Kl) can acquire a small CP even component proportional to e|iir+) through K'^'-K^ 
mixing. This component can then decay into two pions without violating CP. How- 
ever, the CP odd component of \K]^), proportional to \K^), can directly decay into 
two pions if CP is violated in the decay. These two decay modes are illustrated in 

Ki ~ K- H- eA"+ 



indirect CPV e 



direct CPV ^ e' 



TTTT 



TTTT 



The flrst mode of decay is called indirect CP violation and is parametrized by a small 
number e whose relation to e will be given shortly. The second mode is call direct CP 
violation and is parametrized by an even smaller number e'. These parameters are 
defined through 



and 



T[Ks ^ (7r7r)o] 

> (71^)2] ^ T[Ks ^ iTnr)2] 
y/2 \ T[Ks ^ (w)2] ' T[Ks (7r7r)o] 



(t[Kl 



(3.17) 
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A little algebra allows us to relate these CP violating parameters to the measured 
quantities 7700 and rj^ : 



7?oo = e - 2 



11+' 



1 - V2uj 
e' 



(3.18) 
(3.19) 



Moreover, under the assumption that (7r7r)o dominates the sum for in Eq. (3.11), 
the relationship of e to e can be calculated to be 



■ImAo 

' KcAn 



1 + ig-i™-^" 



, ImAo 
ReAo 



(3.20) 



to first non-trivial order in the small CP violating quantities ]^^° 



algebra and some simplifications allow us to express e, e' 



and e. A bit more 
and w in terms of Aq , A2 and 



ImMi2 (de Rafael, 1995) which, in principle, can be calculated using lattice QCD: 



Re^o 



f ImMi2 ImAo 



e<<52-5o) A2 

7= — Iin^ 

^/2 Aq 



(3.21) 
(3.22) 
(3.23) 



where the phase of e is approximately given by (pf^ ~ (2(/)^ + (/'oo)/3 ^ 7''/4, with (p^ 

and (/)oo the phases of rj^ and 7700, respectively. |^In the expression for e', the imaginary 

part of A2/AQ measures the relative reality of A2 and Aq, which is what we need for 
direct CP violation since the latter must arise from the interference between the two 



available decay channels. Note that with the approximations used here, Eq. (3.141 gave 
AMx ~ 2ReMi2, where ReMi2 also could, in principle, be computed on the lattice. 
Experimentally, the various quantities which describe the K^-K^ system and K ^• 



TTTT decays are well measured (Nakamura et al., 2010): 



In many phenomenological studies, 4>e is fixed to 7r/4 and is neglected. However, at the 

levels of accuracy currently reached in the computation of the local contributions to ImMi2, thes e 
approximations are becoming too crude. This has been emphasized in l |Buras and Guadagnoli, 2008[ ), 
where the implications of an estima te of and of t he deviation of ijij from 7r/4 has been inves- 

tigated. Moreover, as explained in (Buras et al, 20101, if , which approximates "^jfj^Jp^i is 

included in Eq. ( |3.22[ l, consistency requires that one also account for long distance contributions to 
ImMi2 . 
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AMk = (3.483 ± 0.006) x lO"^^ [o.2%] 



1 

M 

\e\ 



An 



Re- 



^2 

(20+- 
1 - 



22.4 



(A/ = 1/2 rule) 

- |?7oo|)/3 = (2.228 ± 0.011) x 10"^ 
0oo)/3 = 43.51 ±0.05 [0.1%] 



[0.5%] 



(1.65 ±0.26) X 10" 



[16%] 



(3.24) 

(3.25) 

(3.26) 
(3.27) 

(3.28) 



Using lattice QCD, the weak matrix element relevant for the short-distance, Stan- 
dard Model contribution of |e| has been calculated with a precision of less than 3% (see 



(Lellouch, 2009 Lubicz, 2009 Sachrajda, 2010a Colangelo et al, 2010) for recent re- 



views). We are just beginning to provide phenomenologically relevant information for 
the A/ = 1/2 rule [^0/^2! (|Sachrajda, 2ni0a||Liu, 2010[ [Christ, 2010a|), but Re(e7e) 



is still out of reach for the moment (Sachrajda, 2010a). Moreover, as already men 



tioned and discussed further below, AMk has long distance contributions which make 
its determination on the lattice difficult. However, for that also, progress has been 



made (Christ, 2010&). 



3.2 K^-K^ mixing in the Standard Model 

As we have just seen, K^-K^ mixing arises from the AS" — 2, sd ^ sd FCNC. This 
mixing is responsible for the K^-Ks mass difference, AMk, and indirect CP violation 
in K t:t: decays. 



In the Standard Model, it occurs at one loop through diagrams such as ( Glashow, 



Ihopoulos and Maiani, 1970 Gaillard and Lee, 19746) 









^ 

s ( 


; ^ 


! ^ 

) d 




) < 




d < 


\ U2 < 

■« 1 


) s 









(3.29) 



Setting the external four-momenta to zero, the amplitude associated with this diagram 
is 

(3.30) 

where the u's and w's are the usual particle and antiparticle spinor wavef unctions, 
and we have defined the amplitude with a minus sign for convenience. The W boson 
propagator in Feynman gauge is 



(3.31) 
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and the sum of the up quark propagators, 



^ — niij + it 



(3.32) 



with A(7 = Vt/sFJ^. 

The unitarity of the CKM matrix imphes that X];7=n c t = 0. Thus, we have 
(dropping ie for the moment) 



S{k) = 5] Ac 



U=c,t 



(3.33) 



From this we see the GIM mechanism (Glashow et al., 1970) in action: if to„ 



mt there would be no K'^-K^ mixing. In fact, this process was used to estimate the 
charm quark mass before it was actually discovered (Gaillard and Lee, 19746). 
After performing some Dirac algebra, we find: 



M 



r<2 n/r2 



where we have used Gp — 5|/(4\/2M^). Setting to„ = 0, 



(3.34) 



T, 



4i 



UxU2 



fc2(l-fc2/M2,)2fc2 



(3.35) 



With the spinor wavefunction factors appropriately replaced by quark field opera- 
tors, M can be interpreted as an effective Hamiltonian whose matrix element between 
a K'^ and a K'^ state yields the off diagonal matrix element (3.10) of the mass matrix 
of Eq. (3.7 1. Of course, to obtain the full effective AS' = 2 Hamiltonian, one must 



include the con tributions of all of the diagrams which contribute to the process (jlnami 
and Lim, 198"T| ).p] 



Now, the values of the CKM matrix elements as well as of the masses m„, 
and mt imply that ReMi2 ^ ImAfi2 and that ReMi2 is dominated by the cc term. 
Naively, 



4^2 



^4 



d^k 



ml 

2 



(fc2 + ie)[/c2 — m2 + ie 



O 



1 



O 



(3.36) 



However, a closer look at this loop integral indicates that it is dominated by momenta 
in the range between and mc and we should not forget that all sorts of gluons with 



^Eq. ( 3.34 1 and the amplitudes associated with the other contributing diagrams must be multiphed 
by 1/2 if their spinors factors are replaced by operators to yield the effective /\S = 2 Hamiltonian. 
Indeed, the operator (dL7M'^L)(c'L7'^*L) has twice as many contractions with the four external quark 
states as there really are. These extra contractions correspond to a doubling of the individual box 
diagrams. 



28 A lattice case study: K — )• tttt, CP violation and A/ = 1/2 rule 



or without quark loops can be exchanged between the quarks in the diagram. Since 
these momenta include scales of O(Aqcd) or below, it is clear that corrections are 
out of control and ReMi2 cannot be calculated in this way. 



Said differently, box diagrams, such as the one of Eq. (3.29), cannot be viewed like 



a point interaction, but rather receive long-distance contributions from intermediate 
cc states. In the language of Eq. (3.10), these diagrams contribute to the second term 
through (c?l7^sl)(cl7^ci,) effective operators in 'Has=i- 

The situation is very different for the calculation of the CP violating parameter 
e ~ ImMi2/ReMi2 1^ Indeed, indirect CP violation in X — > vrvr decays comes from the 
interference between the following types of contributions 




and 




So, the relative weak phase between the two types of contributions is arg{(A* )^Mi2} 
and thus, what we really have is 



Im [(A,:)^Mi2] 
Re[(A*)2Mi2] ' 



(3.37) 



where the (A* )^ had been canceled in earlier expressions, using the fact that A„ is real 
in our conventions. 



after Eq. (3.35), we have 



Now, using Eq. (3.34 ) with the replacement of wavefunctions by operators discussed 



J 12 



167r2 



[iXlXt)^Tu + {XlXcfT^c + 2(A:A,)(A:At)r,t) 



where 



v-a\K ) , 



(3.38) 



^For the sake of clarity, we neglect here the small contribution to e from (see Eq. (3.22 i) 



The theory ofK^-K^ mixing 29 



{ds)v^A = dj^il - 75)s 



(3.39) 



and w here th e J'^'^-' are the quartic rephasing invariants defined in Eq. (2.27). As shown 



in Sec. 



2.3.2 



Im/: 



■dcs 



J — — Im/^^^^g, where J is the Jarlskog invariant. Therefore, 



Im[(A:)2Mi2] 



8M^ 



x{K\ds)v^A{ds)v^A\K'') , 



(3.40) 



and J appears as it should in a CP violating quantity. 

Moreover, the integrals in (Ttt — T^t) and {T^c — T^t) are dominated by momenta 
between rric and toj. The same is true of the integrals which appear in the other 
diagrams ( Inami and Lim, 1981[ ) that contribute to this AS* = 2 process. That means 
that the QCD corrections are calculable to the extent that rric can be considered a 
perturbative scale. Thus, under this assumption, and using the experimental value of 
AMjc in lieu of Re [(A* )^Mi2] , we can reliably calculate e with the replacement 



Im < 



Vu^s 
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i.e. by replacing the box diagram with the local, four-quark operator of Eq. (3.40), 
and the appropriate short-distance QCD corrections, omitted here. As mentioned in 
footnote [4j corrections to this approximation have been examined in (Buras et al., 
20T0k. 



3.3 The theory of K^-K^ mixing 

The calculation of M12 of the previous section has actually been performed to leading- 



log (LL) (Vainshtein et al, 1977&| Witten, 1977| Vainshtein et al, 1977 a Oilman and 



Wise, 1983) and next-to-leading-log order (NLL) (Buras et al, 1990) in QCD (for a 
review see (Buchalla et al, 1996)). The resulting NLL, A5 = 2 effective Hamiltonian 
is given by 



167r2' 

xC(^) X {ds)v-A{ds)v-A 



[X^tilttStt + Xl 



(3.41) 



where the running of matrix elements of the four quark operator is canceled by the 
coefficient C(/x), the Sqqi are the Inami-Lim functions (Inami and Lim, 1981) which 
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correspond to the electroweak box contributions in the absence of the strong inter- 
action, and the rjqqi are short-distance QCD corrections. In Eq. (3.41) all Standard 
Model degrees of freedom with masses down to, and including, that of the charm quark 
are integrated out. Thus, in Eq. (3.411, all QCD quantities must be evaluated with 
three active quark flavors. 

It is useful to consider the case of a general number of active quark flavors, Nf, 
and of colors, N^, when discussing the running of this A5 = 2 operator. Indeed, the 
running is the same for other Standard Model, AF = 2 operators, such as AB = 2 
and AC = 2, except that Nf must be chosen appropriately. Moreover, working with 
general Nc allows one to consider the large- A'c limit, to be discussed shortly. With 
as = Q;s/47r, the running of as and of C(/z) is given by: 



din as 
dln/^2 
dlnC 
d In 



-p{as) = -lioas- Pial+0{al) , 
l{as) = loas + 71 -1- 0{al) , 



(3.42) 
(3.43) 



where /3(as) is the QCD /3 function and 7(03) is the anomalous dimension of the 
AS" = 2 operator. It is well known that /3o, Pi, and that LO anomalous dimensions are 
renormalization scheme independent. At two loops, we have (Nakamura et al., 2010) 

nNr~2Nf ^ 34. .0 10 _ „ N?-l 



34 , 10 



and, in the MS-NDR scheme for 71, 



70 = 3- 



Nr 



71 



4iV, 



-21 



57 19 



-Nf 



2Nr, 



It is straightforward to integrate Eq. (3.42): 

/3oln 



1 



(-)'- 


^y 


Kl^oJ 


Po 



as(/io) Po 
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as(Mo) 



{Pl/Po)as{^l) 



(3.44) 



(3.45) 



(3.46) 



l + (/3i//3o)a.(Aio) 



For reasonably small as(/xo) ln(/Lt//xo)^, the coupling at n can be related to the one at 
^0 through: 



asQi) = asit-in) ■( 1 - /3oas(A*o) In 
Pi - Po In 



- asipo)^ In 



0{al 



(3.47) 
(3.48) 



Alternatively, we can define Aqcd as the value of /io at which Os (^0) is infinite, yielding: 
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a, in) 



3i In In 



^^QCD 



(3.49) 



It is also straightforward to integrate Eq. (3.43), 



g(M) 



exp ■ 



(3.50) 



The coefficient C(/i) is only defined up to an integration constant. For consistency 
with Eq. (3.41 1, I consider 



C(/i) = [Anasifi)]-'"'/^" exp. 
Thus, at NLO 

C{^l) = [4^a,(//)]-^°/^° 



""^^''^ da, 
a. 



To 
/3o 



1 + -^a. 
Po 



O(a^) 



20 _ 21 

.-30 /3l, 



(3.51) 



(3.52) 



For the Standard Model AF — 2 operator relevant for K^-K^ mixing, we use the 
anomalous dimension coefficients of Eq. (3.451, with Nj — 3 and, of course, Nc = 3. 
Thus, at NLO, 

307 as{^l)' 



C(M)=a,(M)-^/^ 



1 



162 47r 



Then, in Eq. (3.41 1 the QCD corrections riqq> are of the form 

[l + 0(a,)] 



C(TOc 



(3.53) 



(3.54) 



and where Xq = {mq/Mw) ■ For details, please see (Buchalla et al., 1996). 

To calculate e, we must compute the matrix element {K" \{ds)v-A{ds)v-A{f^)\ K° ) , 
where the kaons are at rest. This is clearly a nonperturbative QCD quantity because 
the typical energies of the quarks within the kaons are on the order of 100 MeV, a 
regime where perturbation theory fails and confinement effects must be taken into 
account. This is where lattice QCD enters. 

For historical reasons, and because it is very convenient in lattice computations, 
we define a normalized matrix element 



_ {K^\isd)v-Aisd)v^Aifi)\K° 



(3.55) 



l{K0\s-f^j5d\0){0\sj^-f5d\KO) 

where we have considered the AS =^ —2 matrix element to conform with convention. 
The benefit of this normalization on the lattice, is that the resulting Bk parameter 

^Here again we choose to neglect the small contribution to e from . Computing it on the 

lattice is a whole other project which is related to the computation of e'. That computation goes 
beyond the presentation I wish to make here. 
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is dimensionless. Therefore it does not suffer from ambiguities due to scale setting, 
something which was particularly bad in old quenched-calculation days. Moreover, 
the numerator and denominator are very similar, and both statistical and systematic 
uncertainties cancel in the ratio. Of course, the convenience of this normalization would 
be limited if the denominator were an unknown, nonperturbative quantity. However, 
the matrix elements in the denominator define the leptonic decay constant of the kaon, 

{Q\-s^^-lM\K\p)) - ilKP^e-'^-"" , (3.56) 

which is well measured experimentally or straightforward to compute on the lattice. 
(In the convention used here, Jk — 156 MeV.) Thus, once the i3-parameter has been 
computed, the normalizing factor is a known quantity and the desired matrix element 
is easily obtained, from 

{k''\ {sd)v.A{sd)v-A{li) \K'') = \flMlBK{ti) , (3.57) 



where M]^ is best taken from experiment. 



Historically, the denominator in Eq. (3.55) was an approximation used to estimate 
the matrix element of 0^^^_2- It is called the vacuum saturation approximation, or 
VSA for short. It is obtained by inserting the vacuum in all possible ways between all 
possible quark-antiquark field pairs formed from the fields of the four-quark operators, 
using Fierz transformations if necessary to bring the fields together. For the case of 
interest here, 

(0|s-,7m(1 - l^)d\x)\K\p)) = -^^J^p^e-'P- , (3.58) 

where a and b are color indices and the dependence on the number of colors Nc = 3 is 
made explicit. Therefore 

(if°|O™_2(M)|if°)y5A-2((i^°|s-,7^(l-75)dl0)(0|s-fc7p(l-75)rf''l^°) 

+ (i^0|s,7'^(l -75)d''|0)(0|5f,7^(l -75)rf'^(a;)|if")) , (3.59) 

where repeated color indices are summed over and where the factor of 2 comes from 
the fact that the two factors of {sd)v-A in 0^^_2 are interchangeable. Plugging 



(3.581 in (3.59) yields 



2^Mi/i . (3.60) 



This is clearly a rather crude approximation, as the LHS of Eq. (3.59) is fi depen- 
dent while the RHS is not. This approximation introduces a renormalization scale 
dependence which is unphysical. 

®Note that the K'^ does not act ually decay leptonically: the do. However, in the isospin limit, 



the decay constant defined in Eq. (3.56 I is equal to the physical decay constant fx of the K^. 
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A more modern approximation to the matrix element is obtained by keeping the 
leading term in a large- A^c expansion. The large- A'c, or 't Hooft limit ('t Hooft, 1974), 
is defined by taking iVc ^ oo while holding a^iVc fixed. By counting the number 
of Q!s and loop factors of Nc in the various contributions to the relevant correlation 
functions (see below), it is straightforward to convince oneself that in the large- iVc 
hmit, i?A'(M) — 3/4. This corresponds to dropping the second term in Eq. (3.59), 
which is clearly suppressed by a factor of l/N^ compared to the first. As in the VSA 
approximation Bk{h) is /J. independent, but here the /i dependence is also absent in 
the short-distance, Wilson coefficient, as can be seen by taking the large- TVc limit in 
Eqs. (3.44)-( 3.51[ ): the large- A*"c approximation is a well-defined and self-consistent 
approximation scheme. 

Before closing this section, it is worth mentioning that one can define a renormal- 
ization scheme and scale independent i?-parameter, B^^, by multiplying Bk{ij,) by 
C(At) of Eq. ([3^ (with iV^ = 3 and iV/ = 3): 



nRGI 



C'iii) X BkM 



(3.61) 



3.4 Computation of bare Bk 

On the lattice, the numerator of Bk is obtained from three-point functions. Quark 
propagators are given by: 

Sq[x,t;ri,ts;U] ^'^D~'^[x,t;Xs,ts;mq;U]ri{xs) , (3.62) 

where D is the lattice Dirac operator associated with the chosen fermion action, to^ the 
quark g's mass, U the gauge field configuration on which the propagator is computed, 
T]{xs) is a three dimensional source which may be a delta function or may have some 
spatial extent and tg is the timeslice at which the source is placed. If only propagators 
from a point source at, for instance, ts — Q and Xg = Q (i.e. rj{xs) = g), are available, 
then we can consider the following three-point function, in Euclidean spacetime of 
course: 

C3{U,tf) = {dl5s{xf,tf)0^^^_Md-^5s{x^,U)) , (3.63) 

X,,Xf 

where the argument of Oas^_2 is its spacetime position, not the renormalization 
scale. For T/2 > -U > 1/AEk and T/2 :^ tf :s> 1/AEk, d'y5s{xi,U) creates a K° 
aX t = ti, this kaon then propagates to t — where O^^_2{0) transforms it into a 
and d"i^s[xf,tf) destroys the resulting at i = t/. /S.Ek is the energy of the 
first excited state in the neutral kaon channel minus Mk- The sums over Xi and Xf 
put the initial and final kaons at rest. Thus, in this limit 

C3[U,tf) — > (0|d75s(0)|i^"(0)) (3.64) 

X (^°(0)|OA5=-2(0)|i^"(0))(i^"(0)|J75s(0)|0) . 

desired mat. elt. 
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This result is obtained by inserting a complete set of hadron states between O|;'g^_2(0) 
and dj5s{xi,ti), between dj5s{xf,tf) and O^^_2{0) and between dj^slxijii) and 
dj5s{xf,tf), keeping the sequence of states which gives rises to the smallest exponential 
suppression and which has the appropriate quantum numbers to give nonvanishing 
matrix elements. If T/2 3> tf,—ti is not realized, then contributions which are not 



significantly exponentially suppressed compared to the one in Eq. (3.64 1 will have to 
be added. 

Similarly, to obtain the matrix elements required to construct the denominator of 



Eq. (3.55), we can consider the following two-point function 



C2Ai) = Y.^d^^<^,t)-s^^^5dm '"^^^ (3.65) 



(p-MKt _ p-MK(T-t)\ 



2M 



K 



(0|d755(0)|X'^(0))(7^"(0)|s757^.rf(0)|0) , (3.66) 

V ' 

denom. mat. elt. 



where I have not assumed here that T/2 ^ t to illustrate the additional contributions 
which arise in that case, and where I have used the properties of the correlation 
function under time reversal. 

Then, Bk is obtained from the ratio 



— > DK{a) , {6.K>i) 



where I have reinstated T/2 ^ tf — ti to get rid of "backward" contributions. In 



Eq. (3.67), the argument a of Bk is there to indicate that this is the value of Bk in 



the lattice regularized scheme and that it still requires renormalization. 

To actually compute C^{ti^tf), we have to take the propagators of Eq. (3.62), 
contract them in the following way and average these contractions over the gauge 
ensemble, i.e. 
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double trace terniocAr? 



s 0^s=-2 d 
t = 



- 75 




(3.68) 



single trace termocAfc 

where Afu is the number of independent gauge configurations. To get the time-reversed 
propagator required to construct the correlation function, one usually uses the 75 
hermiticity which most lattice Dirac operators have, 



Sd[6, 0; f , t; U] = 755d[x, t; 0, 0; U]^^ 



(3.69) 



In Eq. (3.681, one clearly sees how, in the large- A'c limit, only the first contraction 



survives. The second contraction provides a 1/Nc suppressed correction. In practice, 
the two contractions have the same sign, and a cancellation operates in the calculation 
of Bk- 

The method described above is actually a poor way to obtain Bj^ because the 
matrix elements of interest, {K°{6)\Oas=-2{0)\K°{6)) and {0\j5j^d{0)\K"{0)), are 
sampled at only one point on the lattice: at the origin. We would gain a factor of 
roughly (LAf^)^ in statistics if we could sample them over the whole three dimen- 
sional volume of the lattice. Thus, a better way to calculate Bk is to have two zero 
momentum sources for the quark propagators, one at ti and the other at tf. One can 
consider, for instance, wall sources: 



Sq[x,t;W,tw;U] = '^[x,t;Xs,tw;mg;U]'^Ss^,g 



(3.70) 



Then, one constructs the following three-point function 



C3 



d 




— s 



AS=-2{x,t)i 75, 
' d 



(3.71) 



^This is because the longest correlation length in the system is so that regions separated 

by that distance should be reasonably decorrelated. 
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where < and tf > are chosen so that there is a range of t around t = such 
that the correlation function is dominated by the propagation of a zero momentum K'^ 
state between ti and t and the propagation of a zero momentum K'^ state between t 
and t f . The gauge field configurations are usually gauge-fixed on the walls because wall 
sources correspond to meson sources and sinks of the form J2x y'^a.{x,tw)j5S°'iy,tw) 
which are clearly not gauge invariant, except for terms along the "diagonal" x = y. The 
gauge is usually fixed to Coulomb gauge. However, one can also not fix the gauge, the 
result being that the sums over the two quark positions in the sources and sinks reduce 
to a single diagonal sum over the positions after the average over gauge configurations 
is taken. 

Using the wall sources, we also construct the two two-point functions {t\y — ti, tf), 

1 "--^ 

C2At,iw) ^ j^Y^T.^^ ' fJ>7M75. (3.72) 

Then we study the ratio of correlators 

I]pC'2,^(i,t/)C2,p(t, ii) 



as a function of t. For ti <S^t tf, R{t) develops a plateau (see Fig. 3.1) such that 



R{t) BK{a) , (3.74) 

so that Bjiia) is obtained by either averaging R{t) or fitting it to a constant over the 
plateau region. 

3.5 Renormalization of the Standard Model \^S\ = 2 operator 

We are not done, however. As we already stated, simply inserting the lattice operator 
'^AS=-2 i"^ C'3(t) and computing R{t) yields the bare BK{a). This quantity is diver- 
gent in the continuum limit and must be renormalized. And it must be done so in a 
renormalization scheme which matches the one used in the perturbative calculation of 
the short-distance Wilson coefficients. 

It is straightforward to show that the full set of AS — —2, AD = 2 operators of 
dimension d < 6 can be written as: 

Oi = Ofl^_^ = {sd)v-A{sd)v-A (unmix) , (3.75) 
O2.3 = {sd)s-pisd)s^p (unmix, mix) , (3.76) 
O4 5 = {sd)s-p{sd)s+p (unmix, mix) , (3.77) 

where he subscripts {sd)s-p and {sd)s+p are defined in analogy with {sd)v-A in 



Eq. (3.39). In Eqs. (3.75)^(3.77) "unmix" and "mix" refer to the color indices. In the 
"unmix" case, the color indices of the quark-antiquark pairs within parentheses are 
contracted; in the "mix" case, the color index of the quark of one pair is contracted 
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X^/dof=0.79, CL = 0.67 
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Fig. 3.1 Plateau fit to i?(f)/(4L^) on the coarse Nf = 2 + 1, staggered, MILC, 
ami/airih = 0.007/0.05 ensemble. The legend shows the nondegenerate pair of domain- wall, 
valence quark masses making up the kaon in the three-point correlation function. The cor- 
related x^/dof and confidence level of the fit to a constant (red lines) are given in the title. 



Taken from (Aubin et al, 20101 



with the color index of the antiquark of the other pair, and vice versa. Because Oi 
Fierz transforms into itself, the "mix" and "unmix" Oi are the same operator. 

To understand the renormalization patterns of 0^^_2 and the other AS" = —2 
operators, it is useful to consider their transformation properties under various sym- 
metry groups. Because we will only work in massless renormalization schemes, the 
SU (3) L X SU {S)r chiral group is a symmetry which is relevant here. Under the action 
of this group, {sd)v-Aisd)v~A transforms in the (27, 1) representation, i.e. it is a 27 
under SU{3)l and clearly a singlet under SU{3)]i, since it is composed only of left- 
handed fields. It is straightforward to derive the SU{3)l x SU{3)ji representations to 
which the five AS* = —2 operators belong: 



Oi-(27,l) ^^''27^5 1 = 27 



O2.3 
04,5 



(6,6) 6®6 = 27®8®1 



(3.78) 



SU{3 



8® 8 = 27® 10© 10© 8 ( 



1 



In Eq. 



(3.78) I have also worked out the reduction of these representations to the 
R, Eightfold Way, SU{3)v group, for reasons which will be clear 
mix" versus "unmix" feature of these operators has no bearing 



diagonal, V — L 
shortly. Note that the 
on their flavor transformation properties as these features pertain solely to color. 

From this, we see that Oas=-2 the only AS* = —2 — — AD operator of dimension 
6 or less, which transforms as (27,1). Thus, in any regularization which preserves 
SU{3)l X SU{3)fi symmetry (or at least in the valence sector), Oas^_2 renormalizes 
multiplicatively. This includes overlap and domain wall fermions, for sufficiently large 
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fifth dimension. Similarly, Eq. (3.78) indicates that the operator pairs (02,03) and 
(O4, O5) may mix within each pair under renormalization, but not with any of the 
other AS* = —2 ~ — AZ? operators. 

The situation is very different for Wilson fermions. Indeed, the Wilson-Dirac op- 
erator breaks the chiral symmetry of continuum QCD explicitly, down to the vector 
flavor symmetry SU{?>)v- As Eq. (3.781 shows, S'[/(3)y is not sufficient to forbid 



'^AS=-2 from mixing with the four other AS* — —2 — —AD operators, O2,... ,3, under 
renormalization. 

To push things further, we turn to parity. Since parity is preserved by Wilson 
fermions, we can consider separately the renormalization of the parity even and parity 
odd components of the operators. For Bk we are clearly interested in the parity even 
part: 

(i^°|(sd)y_A(sd)y_A|i^°) = {K''\{sd)v{sd)v + (sd)^(sd)^|if") . (3.79) 

So let us begin with the renormalization of the parity even part. 

At this point, it is useful to invoke a discrete symmetry transformation known as 



CPS (Bernard et al., 1985). It consists in performing a CP transformation, followed 



by a switching s ^ d. Note that this vector flavor symmetry is only softly broken by 
the mass terms in the action. Therefore, violations must appear multiplied by factors 
of {vris — md). Under CPS, we have: 



CPS 



'S^pid 



sd sd 



57' d 



-57 d 



(3.80) 
(3.81) 
(3.82) 
(3.83) 



Thus, the parity even components of the AS — —2 operators transform under CPS as 

CPS 



= {sd)visd)v + {sd)Aisd)A 
{sd)s{sd)s + {sd)p{sd)p 
Ot.5 ^ {sd)s{sd)s + {sd)p{sd)p 



Oh 



,+ 
2.3 



CPS ^+ 



(3.84) 
(3.85) 
(3.86) 



All of these operators are CPS eigenstates, with eigenvalue +1. Like SlJ(?))v symmetry, 
CPS does not forbid O^ to mix with ... 3, under renormalization. In fact, there is 
no symmetry which forbids 0\ to mix with the other operators and one finds, in 
practice, that they do mix. We have. 



1=2 



(3.87) 



where Z^{a,ii) is logarithmic ally divergen t in the continuum limit, while the mix- 
ing factors, Zii{a), are finite (Testa, 1998). Using the fact that the values of B^'^^ 
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of Eq. (3.61), obtained in the continuum renormalization scheme and in the lattice 
regularization scheme, must be identical, one easily obtains a formal expression for 
the renormalization constant Z^{a,ji), to all orders in perturbation theory: 

'^-f^) da^ jios) ^ das I'^'jas) 

as Pitts) Jo as fS^'^^ias) 



Zna,,)-c.,{l ^i^^ (3.88) 



'0 



with a^^^ — gQ/(47r)^, where is the bare lattice coupling. Expansions in the bare 
lattice coupling are notoriously poorly behaved, and one can usually improve their 
convergence by expressing them in terms of a renormalized continuum coupling, such 



as QsifJ.) at /i = 1/a for example. Applying this to Eq. (3.881, we find 



Jo asP(as) Jo asP[as) J 

where -t"'*^* is the anomalous dimension obtained by rewriting 7'''* in terms of the 
continuum coupling Og. Obviously this change of variable would make no difference 
to Z^{a,fi) were it to be computed to all orders in perturbation theory. However, at 
finite order this reordering of the expansion may improve the convergence of the series. 



It is interesting to expand Eq. (3.89) to one loop-this could be done for Eq. (3.88) 
instead-to explicitly see the relationship between the coefficients of the expansion and 
the anomalous and beta function coefficients. We find: 

Z+(a, ^) = 1 - ^ (^70 InM' + ^^^-^) + 0{al) . (3.90) 

One clearly sees that the constant 0{as) term knows about the two- loop anomalous 
dimension (it is a subleading log) whereas the leading log term is given by 79. At this 
same order, the mixing coefficients are finite and given by 

zu{a)^z[]'^^+0{a^s) « = 2,...,5, (3.91) 
47r 

where the are constants. 

Now, the SU{3)l x SU{3)]i properties of the operators indicate that the physical 
contributions in {K^^\Oi {a)\K^) are chirally suppressed compared to the nonphysical 
ones: 0{p^) vs 0(1) in counting. Thus, even though the mixing of Oi with 

O2V • • jOs is Us suppressed, this suppression can easily be compensated by an 0(10 — 



20) enhancement from the matrix element ( Babich et al, 2006) . This means that the 



necessary subtractions are delicate and it is preferable to avoid calculating Bk with 
Wilson fermions. 

It is also interesting to study the transformation properties of the parity odd com- 
ponents of Oi, • • • , O5 under CPS. We find 

Of = -2{sd)v{sd)A ^ Of (3.92) 
03-3 = -2{sd)sisd)p ^ -02-3 (3.93) 
0.5 = ■ (3-94) 
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Thus, CPS forbids Oi to mix with O2 3 and O4 5 vanish anyway. We conclude that 
Oi, unhke Oj^, renormahzes muhiphcatively: 



Oi (m) = ^1 (a,M)Oi (a) 



(3.95) 



This turns out to be useful for twisted-mass QCD (tmQCD), which is described in 
detail in (Vladikas, 2010). For instance, twisting (u, d) by an angle a and leaving the 



strange quark s untwisted, we have 

[Ovv+aa]cicb = [[sd)v{sd)v + {sd)A{sd)A 

= COSa [Oyv^+AAltmQCD 



QCD 

i sin a \0 



yA+AVJtmQCD 



(3.96) 



So, picking a = 7r/2 (i.e. maximal twist), the above equation implies the following 
relationship between the renormalized matrix elements: 



{K°\Ovv+AA{lJi)\K°] 



QCD 



VA+AV 



tmQCD 



(3.97) 



Then, since CPS symmetry is only softly broken in tmQCD, {K^\OvA+Av{o)\K'^)ti\iCiCYi 
renormahzes multiplicatively. Thus, by working in tmQCD at maximal twist, one 
can compute Bji{p) performing only multiplicative renormalization (please see Tassos' 
course notes (Vladikas, 2010) for details). 

Though we discussed renormalization mostly in terms of perturbation theory, I 
greatly encourage you to perform this renormalization nonperturbatively, with one of 



the methods explained in Peter (Weisz, 2010) or Tassos' (Vladikas, 20101 course notes. 



You may also want to look into the renormalization procedure put forward in (Diirr 



et al, 2010a Diirr et al, 2010&), which enhances the RI/MOM method of (Martinelli 



et al., 1995) with nonperturbative, continuum running (see also (Arthur and Boyle, 



2010 )). Or if you choose to renormalize perturbatively, you should at least do so to two 



loops to ensure that you have some control over the perturbative series. Of course, one 



may argue that the short distance coefficients in Eq. (3.41) are only known to NLO, 



and that there is no point in doing much better in the lattice to continuum matching. 
Moreover, there are other uncertainties in the relation (3.22) of e to Bk, such as the 
one associated with the error on the determination of |V"c6| or with the neglect of 1/m^ 
and of the ImAo/ReAo corrections (see e.g. (Lellouch, 2009)). However, it is admittedly 
a pity to perform a careful nonperturbative computation of the bare matrix elements 
only to introduce perturbative uncertainties through the matching procedure. 



3.6 Final words on K^-K^ mixing 

Given a preferably nonperturbatively renormalized Bk, it must be matched to the 



scheme used in the computation of the Wilson coefhcients which appear in Eq. (3.41 ), 



and computed for a variety of lattice spacings and quark masses. Then you must use 

^''CPS violating terms are proportional to (r?is — m^), so that Ovv+Aa{o-) can only mix with 
higher dimensional operators which are suppressed by powers of the lattice spacing and the latter can 
only contribute discretization errors. 
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the methods described in, for instance, (Diirr et aL, 2008 Lellouch, 2009[[Durr et al, 



2010a[ purr et al., 20106[) and /or in Pilar ( [Hernandez, 2010[ ), Peter ( |Weisz, 2010[ ) and 
Maarten's ( Golterman, 2010| ) course notes to extrapolate (preferably interpolate) to 
the physical values of mud, and extrapolate to the continuum limit. Finally, you 
must perform a complete systematic error analysis, such as the ones performed in 



(Diirr et aL, 2008 


Diirr et al., 2010c 


Diirr et al., 2010a 


Diirr et al, 20106 


). For a 


recent, comprehensive review of lattice calculations of B^, see for instance ( 


Colangelo 



et al., 2010). 



Before concluding this discussion of K^-K^ mixing, it is worth mentioning that 
lattice QCD can also provide information that is relevant for this process beyond the 
Standard Model. Quite generically, when one adds heavy degrees of freedom to those 
of the Standard Model, such as in supersymmetric extensions, one finds that the full 
set of /S.S = —2 operators of Eqs. (3.75 1-( 3.77) contribute to the low-energy effective 
Hamiltonian. Of course, the lattice can also be used to compute the matrix elements 
of the four additional operators, between and states. This has been studied in 
the quenched approximation in (Donini et al., 1999 Babich et al., 2006). Ref. ( )Babich| 



et al., 2006) finds ratios of non Standard Model to Standard Model matrix elements 



which are roughly twice as large as those in (Donini et al., 1999). As explained in 
(Babich et al., 2006), this is most probably due to discretization errors present in 
(Donini et al., 1999). This picture appears to be confirmed by the preliminary Nf ~ 2 
results of (Dimopoulos et al., 2010). 



3.7 Phenomenology of the Al =1/2 rule 



The goal here is to compute nonleptonic weak decay amplitudes, such as those for 
K — > TTTT, directly in the Standard Model, without any extraneous model assumptions, 
nor potentially uncontrolled approximations such as LO, SU{3) xPT. This is critical 
for showing that QCD is indeed responsible for surprising phenomena such as the Al = 
1/2 rule, or deciding whether New Physics is hidden in the experimental measurement 
of direct CP violation in — tttt (i.e. of e'). In the following 1 will concentrate on the 
Al =1/2 rule as the calculation of e' is significantly more difficult. For the latter, the 



renormalization of the relevant matrix elements is more complicated (Dawson et aL, 



1998). There are more matrix elements involved and there appear to be important 



cancellations between them (see e.g. (Buchalla et al., 1996)). Moreover, final-state 
interactions, such as the ones discussed in (Pallante and Pich, 2001 Buras et al., 2000 ), 



seem to play an important role. Indeed, recent attempts to exhibit the Al = 1/2 rule 
and to determine Re(e'/e) using soft pion theorems appear to fail due to the large 
chiral corrections required to translate the K ^ n and K ^ amplitudes computed 
on the lattice into physical K tttt amplitudes (Li and Christ, 2008 Sachrajda, 



2010a Christ, 2010a 



Experimentally, the partial widths measured in the different K ^ tttt decay chan- 
nels, together with the corresponding isospin changes between the final two-pion and 
initial kaon states, are (Nakamura et al., 2010): 



42 A lattice case study: K — > tttt, CP violation and A/ = 1/2 rule 



V 



+_ = r(ifs ^ TT+TT") = 5.08 X lO-i^MeV 



2' 2 



Too = r(if5 ^ 7r°7r°) = 2.26 x lO^^^MeV 
r+n = T{K+ tt+tt") = 1.10 x lO^^MeV 



A/ 



A/ 



1 3 

2' 2 
3 



Using these results, we find 



^= 463. and ^ 

J- +0 J- +0 



205. 



(3.98) 
(3.99) 
(3.100) 

(3.101) 



whereas /r+o should be ^ 1 in the electroweak theory and in the absence of the 

strong interaction. Now, in terms of the amplitudes, the rates are 



r+- = 

Too — 



2V2Re (^Aq 



r+o = J7l^2p 



(3.102) 
(3.103) 
(3.104) 



with 7 = y/M]^ - 4M2/(167rM|-). Considering r+_ +roo and T+o, and taking 
134.8. MeV and Mk = 494.2 MeV (i.e. isospin hmit values), we find 



\Ao\ = 4.66 X 10"'' MeV and |^2| = 2.08 x 10"^ MeV 



(3.105) 
(3.106) 



whereas the combined chiral and large- TV^ limit predicts (Lellouch, 2001 1 ! It is the 



huge enhancement of over 400 in the rate or 20 in the amplitude which is known as the 
A/ =1/2 enhancement. It was termed a "rule" because this enhancement of A/ =1/2 
over A/ = 3/2 transitions is also observed in other decays, such as A — >■ Ntt, S — > Ntt 
and S — > Avr. 



3.8 The A/ = 1/2 rule in the Standard Model 

At leading order in the weak and strong interactions, A5 = — 1, AZ3 = 1 transitions 
occur through the tree level diagram 




u, c, t 



q' = u,c,t 
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As usual, we integrate out the heavy W boson and t quark. At leading order in 
QCD, this yields the following four-quark, local vertex 




q = u,c 



q = u,c 



with 



Ql = {sq)v-A{q d)v- 



A ■ 



(3.107) 



Now, if we include as corrections ( [Gilman and Wise, 1979 ), such as 
s d 




we generate a new operator: 



{sd)v-A{q'q)v^A 

These corrections also lead to the penguin diagram 

ti w 



(3.108) 




with q = u,d, s,c,b. However, one finds that they yield tiny contributions to the 
CP conserving parts of i^T —j' tttt decays. Thus, we neglect these here. Moreover, we 
are only interested in external states with u, d and/or s quarks. So we do not need 
operators with a single charm leg. To simplify the operator structure, we can also use 
the unitarity of the CKM matrix: 



VcdV* 



-V^dV* - VuV* 



(3.109) 



where the second term on the RHS of this equation can be neglected because it is mul- 
tiply Cabibbo suppressed by a factor of A* ~ 0.003 compared to the first term. Then, 
for CP conserving A5' = — 1, AZ? = 1 transitions, we have the effective Hamiltonian 



n 



AS=-1 
CPC 



V2 



(3.110) 



with 



0± = [{su)v-Aiud)v-A± {sd)v-A{uu)v-A] - [u ^ c] , (3.111) 

where the second term, in which u is replaced by c with an overall minus sign, imple- 
ments the GIM suppression mechanism (Glashow et al, 19701. It is straightforward 
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to show that 0+ is in the (84, 1) representation of the SU{4)l x 5'J7(4)fl chiral group, 
appropriate for classifying operators composed of u, d, s and c quarks for renormal- 
ization in massless schemes: it is a completely symmetric and traceless tensor in two 
fundamental and two conjugate SU (4) l indices. 0_ , on the other hand, is in the (20, 1) 
representation of this group (it is a completely antisymmetric and traceless tensor in 
two fundamental and two conjugate SU{A)l indices). Thus, 0+ and 0_ do not mix 



1981): 



under renormalization. In Eq. (3.110), the short distance Wilson coefficients are given, 
at 0{as), by ( Gaillard and Lee, 1974a[ Altarelh and Maiani, 1974 Altarelli et al, 



An 



7i"^ In 



with, in the MS-NDR scheme: 



(0) 

7± 



= ±3 



and 5+ = ± 



11 iVc T 1 

2 Nr_ 



(3.112) 



(3.113) 



One often pushes the short distance calculation further and integrates out the charm 
quark. But the GIM mechanism is very useful for reducing the divergences of relevant 
weak matrix element on the lattice. (Unfortunately, we will not have the time to cover 



this here, but I refer the interested reader to (Dawson et al., 1998), for instance, for a 
discussion of this and related issues.) 



show that 0_ 



A straightforward analysis of the isospin structure of the operators of Eq. (|3.111 ) 

3 [11] o 



is pure 7 = |, while 0+ has both 7 = i, | components. ^Une may 
wonder then, whether the Al =1/2 enhancement comes from the running of the 
Wilson coefficients C±(/i) from the scale /i ^ Mw, where the ratio C-{^) / C+{ijl) 
is 1 plus small corrections of order as{Mw) ~ 0.1, down to a scale ^ ^ 2 GeV. 
This would mean that the A/ =1/2 rule could be understood as a short-distance 
enhancement. At l eadin g-log order (LL) ( Gaillard and Lee, 1974a[ Altarelli and Maiani, 
1974| ), using Eq. (Esol, we find C_ (2 GeV)/C+(2 GeV) = (as(2 GeV)/a,(m6))'**/''^ 
{as{mb) / as{Mw))^'^^ (C_(Afw)/C+(Mvy)) ^ 2. So there is some short-distance 
enhancement, but not enough by a long shot to explain the A/ =1/2 rule. In turn, 
this means that most of the enhancement in 



A2 



((7r7r)o|C+0+-fC_0_|i^°) 
((7r7r)2|C+O+|if0) 



(3.114) 



must come from long distances. 

One may also wonder what role the charm quark plays in these decays. In partic- 
ular, one might consider an imaginary world in which the GIM mechanism is exact, 
i.e. rric = to„. Does the A/ =1/2 enhancement persist in that limit? One way to 



^^In terms of SU{3)l x 5(7(3)^ representations, 0_ is pure (8, 1) while 0+ contains both (8, 1) 
and (27, 1) representations. 

^■^If you have been reading these notes carefully, you will be quick to point out that this statement 
has its limitations. Indeed, beyond LLO, the running of the Wilson coefficients is scheme dependent. 
However, it is difficult to justify not going beyond that order at scales /i ~ 2 GeV. So the statement 
should be be understood as applying to commonly used schemes. 
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address this problem is to work in the S'C/(4) chiral hmit and determine the LECs 
corresponding to 0_ and 0+, i.e. the (20, 1) and (84, 1) couphngs (see comments after 



Eq. (3.1111) (Giusti et al, 20071. Numerically, in the quenched approximation, it is 
found that there is an enhancement of IA0/A2I over naive expectation (e.g. large- iVc), 
but that this enhancement is roughly a factor of four too small. 

An interesting way to understand this SU{A) chiral limit enhancement is to con- 
sider the three-point function contractions required to determine the matrix elements 
(7r+|0±|/f+). It is straightforward to see that in the case of 0+, with the appropri- 



ate flavor replacements, the contractions are the same as those of Eq. (3.681, up to 



the overall factor of 2. On the other hand, for the pure AI = 1/2 operator 0_, the 



contractions are those Eq. (3.68), but with a plus sign between the double and single 



trace terms. Thus, in the large-A^c limit, the two m atrix e lements coincide and we 



find the value of IA0/A2I = discussed after Eq. (3.1061. However, for finite iV^, 
the single trace term contributes, and does so with opposite signs to {n'^\0+\K'^) 
and to {Tr^\0-\K~^) . Since both trace contractions are positive, the larger the single 
trace term, the larger is {Tr~^\0-\K~^) and the smaller is {tt'^\0+\K'^). This creates 
a A J =1/2 enhancement in which both the numerator |^o| is enhanced and the de- 
nominator \A2\ is depressed. The argument also implies an anticorrelation between 
Bk, the i3-parameter of K'^-K^ mixing, and the AI = 1/2 enhancement of IA0/A2I. 
Indeed, in the chiral limit, the smaller Bk is compare d to its large- A^^ value of 3 /4, 



the larger |Ao/^2| is compared to \/2, as first noted in (Pich and de Rafael, 1996). 



Given the argument which we made earlier, one would think that LQCD is well 
suited to study K — >■ tttt decays. However, there are many conceptual problems one 
encounters when trying to study these decays on the lattice. Some of these are: 

• The renormalization of the AS = — 1 effective Hamiltonian is difficult on the 



lattice, even more so for the CP violating part (see e.g. (Dawson et al, 1998)); 

• lattices are only a few fermi in size and the final-state hadrons cannot be separated 
into isolated, asymptotic states; 

• only approximately evaluated Euclidean correlation functions are available 

There are also technical challenges. For instance, the study of these decays requires 
the calculation of 4-point functions. Moreover, power divergences must be subtracted, 
if the charm quark is integrated out in the case of CP conserving decays, and once 
the W and t are integrated out in the CP violating case. Both these points make the 
study of if — > TTTT decays very demanding numerically. 



3.9 Euclidean correlation functions and the Maiani-Testa theorem 

For well over a decade, it was believed that if —> tttt amplitudes could not be studied 
directly on the lattice. Indeed, these amplitudes have both real and imaginary strong- 
interaction contributions while, in the Euclidean, correlation functions are purely real 
(or imaginary). Thus, it was difficult to see how such amplitudes could be extracted 
from a lattice calculation, necessarily performed in the Euclidean. 



Of course, the Osterwalder-Schrader theorem (Osterwalder and Schrader, 1973 



Osterwalder and Schrader, 1975 ) guarantees that Euclidean correlation functions can 
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be continued to Minkowski space, at least in principle. However, in practice such 
analytical continuations are essentially impossible with approximate Euclidean results. 
This led Maiani and Testa, in 1990, to investigate what can be extracted from 



Euclidean correlation functions without analytical continuation (Maiani and Testa, 



1990). They considered the following type of center-of-mass frame Euclidean correla- 
tion function, 

{TT(pMX-pM)'H^i='^K\Q,U)) , (3.115) 

and asked the question: what information does this correlation function contain re- 
garding physical K — > tttt decay amplitudes in the usual lattice, asymptotic limit, i.e. 
ii,i2 > l/M^, -U > 1/Mk? 

To "simplify" the problem and disentangle Euclidean from other possible lattice 
effects, they chose to work in a large, quasi-infinite volume. This apparently innocuous 
assumption has rather important consequences. For one, in infinite volume, the tttt 
spectrum is continuous. This means that in the limit —ti 3> l/Mn and ii, t2 3> 1/M^, 
only the ground state contribution can be picked out: there is no known numerical 
technique to isolate an excited state in a continuous spectrum. In turn, this implies that 
one can only extract information about the matrix element (7r(0)7r(0)|'H^^~~^|i4'(0)), 
in which all mesons are at rest: the physical decay is not directly accessible on the 
lattice. 

This statement became known as the "Maiani- Testa theorem." It was a formal- 
ization of the general belief that K ^ tttt decays could not be studied directly on 



the lattice but rather that approximations (Bernard et al., 1985) or models (Ciuchini 



et ai, 19961 were needed to obtain information about physical K — >■ tttt decays from 
lattice calculations. 

However, as is often the case with "no-go" theorems, the solution is found by 
questioning the underlying, apparently innocent assumptions. Here, it was the infinite 
volume assumption that brought in all of the difficulties while, for simulations per- 
formed in boxes with sides of a few fermi at most, it does not even approximately 
hold. 



3.10 Two-pion states in finite volume 

In a finite box with sides L,p^two pions cannot be isolated into noninteracting asymp- 
totic states. Rather, the tttt eigenstates are the result of a stationary scattering pro- 
cess. In addition, boundary conditions generically imply that the particles' momenta 
are quantized. For periodic boundary conditions, they come in discrete multiples of 
27r/L: k = z(27r/L), with z G Z"^. In turn, this means that the spectrum is discrete 
and the splitting is actually rather large for box sides of a few fermi. In such boxes, 
the typical spacing between momenta is Ap — 2tt/L = 1.2 GeV/L[fm]. This is clearly 
quite different from the continuous spectrum found in infinite volume. 

In the free theory, in the center-of-mass frame and in the Af (i.e. cubic spin-0) 
sector, the energy of the n-th excited state is given by: 



^^L must be large enough so that the range of the interaction between the two pions is contained 
within the box. 
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Fig. 3.2 The two-pion spectrum in a box of volume as a function of L in the I = J — 
channel, under the four-pion threshold W/Mt^ = 4. The free spectrum is depicted by the 
green dashed curves. The solid blue curves show the interacting spectrum as obtained from 



Eq. (3.1211, using the one-loop I — J — phase shift from (Gasser and Meissner, 1991 



Knecht et al, 19951 



(3.116) 



for n < 6-there is no integer three-vector with norm squared, 7, nor 8 for that matter. 
This spectrum is shown as the dashed lines in Fig. |3.2| 

In the presence of interactions, the energy Wn of the fully interacting state was 



worked out by Martin to all orders in relativistic quantum field theory (Liischer, 1986 



Liischer, 1991 ) for two-pion energies below the four-pion threshold, up to corrections 



which fall of exponentially with the box size (i.e. up to finite-volume, vacuum polar- 
ization effects). 

The ground state n = requires special treatment. In the isospin / = and spin 
J = channel, we have (Liischer, 1986): 



0{L- 



where 



ci = -2.837297, C2 = 6.375183 , 
where the S'-wave scattering length in the appropriate channel is a/ with 

aj ~ iim — ; — , 
fc^o k 



(3.117) 
(3.118) 

(3.119) 



and where k is the momentum of the pions in the center-of-mass frame. 
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For excited states, the results are given in terms of the scattering phase, (5/, of the 
relevant isospin channel. The appearance of scattering phases should not be surprising. 
We know from scattering theory that, under reasonable conditions, potentials can be 



reconstructed from these phases. Martin finds that (Liischer, 19911 

Wr, = 2^ Ml + kl n= 1,2,3,... , 
where fc„ is a solution of the quantization equation 

riTT - 5i{kn) = 0(g„) , 

with Qn = knL/{2iT) and with 



tan 0(g) = 



(3.120) 
(3.121) 

(3.122) 



(j){q) is defined for g > 0. It is such that 0(0) — and it depends continuously on q. It 
is given in terms of the zeta function of the Laplacian 



Att ^ 



(3.123) 



for Re s > 3/2 and by analytic continuation elsewhere. A useful integral representation 
for evaluating ZQQ(l;q'^) numerically is given in Sec.|4] 
Solving Eq. (3.1211, one generically finds 



27ry 



O 



L3 



(3.124) 



where the equation actually gives the whole tower of 1/L corrections once the scat- 
tering phase is specified. As already stated, the equation holds for n < 6 and there 
is no integer three-vector whose squared norm is 7 or 8. Then, beginning at 9, there 
are integer vectors, z, which are not related by cubic rotations but which have the 
same squared norm, e.g. zl = (2, 2, 1) and Z2 = (3, 0, 0). In a treatment where only the 
0(3) spin-0 component of the representation is taken into account the states 
associated with such three-vectors in the free case do not feel the interaction in the in- 
teracting case either, and have k = \z\{2tt/L). This is because they combine into states 
with 0(3) spin-4 and/or higher spins. Such states will be ignored in the following and 
we will consider only n < 6, which is certainly not a limitation in practice. 

The full solution (taking for instance / = J = 0) is shown in Fig. 3.2 The first 
remark which can be made is that the two-pion spectrum on lattices which can be 
considered in the foreseeable future is far from being continuous. The second is that 
distortions due to interactions are quite small: the volume suppression of the correc- 
tions is effective when L > 3 fm. Finally, it is clear that by studying the energies Wn{L) 

^^Because the cube is not invariant under generic rotations, the irreducible representations of the 
cubic group are resolved into many irreducible representations of 0(3). For the cubic representa- 
tion, the relevant spin representations are spin 0, 4,. . . . 
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as a function of box size L, we can turn Eq. ( 3.121[ ) around and reconstruct, at least 



for a few discrete momenta, the scattering phase Sj(k). 

Now, suppose that there is a single resonance R in this tttt channel, with mass 
Mil < 4Af^, i.e. under the inelastic threshold, and width F^j. To explain how Martin's 
equation works in that case, it is useful to turn to quantum mechanics. A quantum 
mechanics approach is actually justified because corrections suppressed exponentially 
in L are neglected in the derivation of the quantization formula. Such suppression 
factors generically correspond to tunneling phenomena. Here they are associated with 
features of relativistic quantum field theories which are absent in quantum mechanics: 
the exchange of a virtual particle around the box. 

We begin by decomposing the total Hamiltonian H of the two-pion system into a 
free part Hq and an interaction Hint- 

H = Ho + H,nt- (3.125) 

Then we consider the n-th free tttt state (n < 6), \tio), in our chosen channel (here 
isospin / = or 2 and J'' = 0~). It is such that {no\no) — 1 and 

iloino) - W^Itio) . (3.126) 

To understand what this energy becomes in the presence of interactions and of 
the resonance, it is useful to consider a perturbative expansion in the interaction Hint, 
though the final result is accurate to all orders. Denoting the resulting energy Wn, and 
\n), the corresponding fully-interacting eigenstate, we have, to second order in Hint- 

Wn = {n\H\n) 

^ wr - 

Here a runs over 2, 4, ... pion states, as well as any other state which appears in the 
given channel, and Wa is the corresponding energy. Factors of the form (wl°'* — Wa) 
also appear in the higher-order terms of the perturbative expansion represented by 
the ellipsis. 

As long as L is such that Wa — Mr is far from the free two-pion energy, wii'\ 
all of the terms in the perturbative series are regular and can therefore be resummed. 
Now, if the resonance is narrow, the coupling of the resonance to the tiq, tttt state will 
be small. In turn, this means that the matrix element {'K'K\Hint\R) is small compared 
to the mass and typical energies of the system. Moreover, the leading correction 
which the resonance brings to the interacting two-pion energy, VF„, appears at second 
order in the expansion. Thus, it is of order |(7r7r|_ffi„t|i?)p and is therefore small. 

When L — such that Wn\Lji) = Mr, the effect of the resonance is radically 
different. Its second and higher order contributions to Wn blow up. In such a situation, 

^^We assume here that the narrowness of the resonance is not only due to phase-space suppression. 
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Fig. 3.3 The two-pion spectrum in a box of volume L'^ as a function of L in the 
I = J = Q channel with an additional, fictitious resonance of mass Mr — Mk^ and 
width Vr = 10^^ X (r+_ + Too - (4/3)r+o) ~ 7.3 MeV with r+_, Too, and T+o given in 



Eqs. ( 3.102 1-(3.104 1. The dashed green curves represent the two-pion spectrum interacting 



through (5o (they correspond to the solid blue curves in Fig. 3.2 1. The solid blue curves show 



the interacting spectrum as obtained from Eq. ( |3.121[ ), with the contribution of the resonance 
to the phase shift. At the points at which the resonance (i.e. the horizontal green dashed curve 
at Mks/Mtt ~ 3.57) crosses the excited I = J = states, one clearly sees a level repulsion 
effect. The effect is rather small and limited to a smallish region around the crossing point 
because the resonance is narrow: Tr/Mr ~ 1.5%. 

we have to resort to degenerate perturbation theory and first diagonalize H in the two- 
state subspace {|no)j This means diagonaUzing the 2x2 matrix 



{no\H\no) {n„\H\R) 
{R\H\n„) {R\H\R) 



( Mb. M„ 

\M*^ Mr 



(3.128) 



where M„ = {'no\H\R) is the transition amplitude between the resonance and the 
two-pion state |no). A straightforward diagonalization yields 



W^=Mr±\M^\ 



(3.129) 



thereby lifting the degeneracy and giving rise to a typical level repulsion phenomenon. 



Thus, in solving Martin's formula (3.121 1, we would find a dependence of the two-pion 



energy as a function of L, which looks like what is depicted in Fig. |3.3| 



3.11 — )• TTTT in finite volume 



What Martin and I realized over ten years ago is that to study K — )■ tttt decays in 
finite volume, we could treat the kaon as an infinitesimally narrow resonance in the 
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weak interaction contribution to the scattering of the two pions ( Lellouch and Liischer, 



2001). In that case, one considers the free Hamiltonian to be the QCD Hamiltonian, 



I.e. 

Ho = i^QCD , (3.130) 

and the perturbation, Hinn to be the effective weak Hamihonian relevant for K — > tttt 
decays, i.e. 



H,nt = Hw= d'x Hw{x) . (3.131) 

J xo=0 

Then, since the amplitudes for K — > tttt decays, T{K — > tttt), are computed at 0{Gf), 
we can perform all of our computations to that order. 

Following what was done in the preceding section for the resonance, we tune the 
size of the box to L = Lk, such that for some level n, 

WniLK) = Mk . (3.132) 

In that case, the corresponding pion momentum is the momentum, k,^, which the pions 
would have in the physical kaon decay, i.e. 



K{LK) = K^\j^ -Ml . (3.133) 

Then, the transition matrix element in the finite volume V = L^, 

Ml^= v{{TTTi)in\Hw\K)v , (3.134) 

is an energy conserving matrix element. Since finite-volume corrections to a single, 
stable particle state are exponentially small in L and since we neglect such corrections 
here, \K)v is identical to jif), up to a purely kinematic normalization factor. 

However, we are not interested in the finite-volume matrix element of Eq. (3.134). 
What we want is the infinite-volume transition amplitude, 

Ti EE {{'K'K)in,out\nw\K) , (3.135) 



where the corresponding Aj of Eq. (3.3) can be made real in the CP conserving case. 

Because it is important here, let us pause to say a few words about the normal- 
ization of states used. In finite volume we use the usual quantum mechanical nor- 
malization of states to unity. Thus, for a spinless particle of mass m and momentum 

v{p\p)v = • (3.136) 
In infinite volume, it is the standard relativistic normalization of states, i.e. 

(pb')-2pO (2^)3^(3) (p-^p-) , (3.137) 



with p*^ — \J vn? + p^, which is implemented. 

To obtain the relationship between these two amplitudes, we compute the shift in 
energy brought about by the presence of the weak interactions in two different ways, 
and require the two results to agree. 
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Fig. 3.4 Kaon contribution to tlie elastic weak scattering of two pions in tiie S'-cliannel. The 
shaded vertices represent the K — > tttt transition amplitudes at 0{Gf)- 

We begin by repeating the degenerate perturbation theory of the preceding section. 
We obtain 



W;^^MK±\Mi\ 



(3.138) 



where is clearly 0{Gf)- Then, we turn to Martin's finite-volume quantization 
formula (3.121). In the presence of the weak interaction, the phase shift receives a 



weak contribution, 6w- Thus, in that formula, we have to perform the replacement 



Sj + 5w . 



(3.139) 



We are interested in this phase shift at the values of momenta k = , correspond- 
ing to the perturbed energies of Eq. (|3.138|): 



Wn\Mi\ 

4:kr, 



OiG%) 



(3.140) 



Indeed, we know that those values must come out of the quantization formula because 
the two methods of determining the energy shifts must give the same result. Because 
are "infinitesimally" close to Mk (in our LO counting in G_f), 6w{kn) is dom- 
inated by the s-channel kaon exchange depicted in Fig. |3.4[ Any other contribution 
will be at least 0{G^). In the s-channel, however, the factor of Gj^ coming from the 
two K-mr vertices is compensated by a propagator enhancement, due to the fact that 
we are sitting only 0{Gf) away from the peak of the resonance. Indeed, the scattering 
amplitude corresponding to Fig. |3.4| is 



5[Fig.|3;4] = 




2H^r|M,{| 



O(G^) , 



0{Gy) 



(3.141) 



(3.142) 



where we have used the fact that the vertices in Fig. |3.4| are the on-shell transition 
amplitudes up to higher order corrections in Gp. To translate this amplitude into a 
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scattering phase, we use the partial wave decomposition of a tttt scattering ampUtude 
S in the center of mass frame: E3 



S = 1(S^WY,{21 + l)P,(cos e) e^''^^^ sin5[(fc) 



(3.143) 



where W = 2y^ + fc^ is the energy of the pion pair. Since the ampUtude of 
Eq. (3.1411 has no angular dependence, it only leads to a zero angular momentum 



phase shift, which is the weak phase shift of interest. Thus, we find 



327r(W^i°))2|M^| 



0{Gl) . 



(3.144) 



We can now include this contribution into the total phase shift 5i and write down the 
resulting quantization equation: 



riTT — 5i{kn ± Afc) ± 



kt\Ai? 



32^(M^i°^)2|M„^| 



± Ag) + 0{Gl) , 



(3.145) 



with Aq = LAfc/(27r). Expanding this equation to 0{Gf), we finally find the relation- 
ship between the desired, infinite- volum e amplitude Ai and the fin ite-volume matrix 
elements M^, computed on the lattice (Lellouch and Liischer, 2001 ): 



d 



d 



were we have used the definition 



fc=fc„ 



Mi = v{{7rTr)in\nw{0)\K)v = M^/L^ 



(3.146) 



(3.147) 



In our derivation w e have t uned th e size of the bo x to Lj f such that Wn^^ — Mx and 



thus, fc„ = k^ (Eqs. ( |3.132l )-( |3.133[ )). However, Eq. ( |3.146[ ) is also vahd for ^ M^, 
as can be seen in Sec. 3.12|and as was derived using the fact that the matching factor is 



related to the density of interacting two-pion states in finite volume (Lin et at, 2001 ). 



An interesting discussion of this and other ways of looking at this formula is given 



in ( Testa, 2005, ). In [Kim et ai, 2005 ), Eq. (3.146) was further generalized to moving 
frames, i.e. frames in which the center of mass has a nonvanishing momentum. And in 



(Kim and Sachrajda, 2010), it is shown how partially-twisted boundary conditions can 



be used to obtain the phase shift (52 (fc) and its derivative in the isospin-2 channel. With 
twisted boundary conditions, one allows some of the quark flavors to be periodic only 
up to a phase. This phase forces the flavors concerned to carry a momentum which is 
proportional to the phase. The boundary conditions are called partially-twisted when 
it is only the valence flavors which are given a twist. 



^^Note that S^'s subscript in Eq. (3.143 i corresponds to angular momentum I, while the superscript 
is the isospin /. This notation will be used in this equation only. Elsewhere, S's subscript will be the 
isospin, except for Syy where W stands for weak. 
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The proportionality factor (3.146) is to a large extent kinematic, as it accounts 
for the difference in normalization of states in finite and infinite volumes, given in 
Eqs. (3.136 )-(3. 1371. This can easily be seen in the absence of interactions. To reach 

the cube must have sides 



the n two-pion with energy Wn 



(0) 



and pion momentum fc„ 
27r ^ 



(0) 



Z,(0) 



Then, Eq. (3.1461 assumes the form 



/|2 



where 



number of z G Z 3 



(3.148) 



(3.149) 



(3.150) 



The proportionality constant is precisely the relative normalization of free kaon and 
two-pion states in finite and infinite volume projected onto the and spin-0 sectors, 
respectively (see also Sec. 3.12). The constant is the product of the ratio of squared 
norms of the kaon state, 2MkL^, and of the two-pion state, 2{WnL^)'^ , times the 
square of the factor relating Hw and 'Hvy(O), i.e. times l/f„ since the finite- 

volume, state is obtained by summing over the pion momentum directions. 

3.12 K — )■ TTTT in finite volume: a simple relativistic quantum field 
theory example 

To understand how the finite-volume effects predicted by Eqs. (3.121 1-(3. 146) show 



up in correlation functions similar to those one would use in numerical simulations, 
it is useful to consider them in the context of a relativistic field theory in which 
all quantities of interest can be computed analytically. Because the form of the finite- 



volume formulae (3.121 ) and ( 3.146 ) does not depend on the details of the dynamics, we 



choose to work in a world in which this dynamics is simplified, so as not to obscure the 
discussion of finite-volume effects with superfluous technical details. The calculations 



below were summarized in (Lellouch and Liischer, 2001). 



3.12.1 Specification of the model 

We consider a theory of a single, neutral, spinless pion field Tr{x), of mass M^. In the 
notation of Sec. |3.lH the unperturbed Hamiltonian density is 

Ho=Hki„+Hs, (3.151) 

where T^kin is the usual kinetic Hamiltonian of a scalar field and the "strong" interac- 
tion between the pions is given by 

A . 



(3.152) 



We assume here that the theory is perturbative in A and we will work to first nontrivial 
order in A, i.e. 0(A). Because the lattice calculations are performed in Euclidean 



^'^The factor of 2 is required because the two pions in the final state are identical particles in the 
isospin limit which we consider here. 
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spacetime, we rotate this theory into the Euchdean and consider Euchdean correlation 
functions. Moreover, our calculations will be performed in a three-volume with 
periodic boundary conditions, but the time direction will be considered of infinite 
extent. 

To make the perturbation theory completely well-defined, we introduce a Pauli- 
Villars cutoff A. At tree level the Euclidean pion propagator is then given by 

SA^) = I e^'-i^i^MO)) = ^ - ^ . (3.153) 

The cutoff should be large enough so that ghost particles cannot be produced at 
energies below the four-pion threshold, but in view of the universality of Eq. (3.1211 
and (3.146) there is no need to take A to infinity at the end of the calculation. 

Since we are going to be computing correlation functions in the time-momentum 
representation, it is useful to have the pion propagator in this same representation. 
We have 

S^{t;k)= [ e-*^~(^(i,f)^(0)) 



-1^ e-^^ N - e-^^ 1*1 , (3.154) 



where Ek ^ JP + and £k = VP+A^. 

As far as the kaon and its decays into two pions are concerned, the least complicated 
possibility is to describe it by a free hcrmitian field K{x) with mass Mk and to take 

nw = I^TT^ (3.155) 

as a weak Hamiltonian density. We will only work here to leading order in the weak 
coupling, also. 

3.12.2 Determination of the phase shift 



Let us first determine the phase shift, S{k), in the model of Eq. (3.152). This is an 
infinite- volume, Minkowski space calculation, though at the order at which we work 
this fact makes very little difference. The partial wave decomposition of the invariant, 
scattering amplitude S, in the center-of-mass frame, is given in Eq. (3.143). At 0{X), 
S — —X, and therefore 

^ik) = -~+OiX^), (3.156) 
where W = 2Ek is the free two-pion energy. 



3.12.3 Two-pion energies using Eq. (3.121) 



In the absence of interactions, i.e. for 5{k) = 0, the solutions of Eq. (3.121), for 
n = 1, 2, . . . , 6, are the free, finite- volume momentum magnitudes, kn^ = y'n{2Tr / L) 

^^Here and in the following, quantities with the superscript (0) are computed at 0(A"), while 
thos e with out a superscript are the same quantities in the presence of the "strong" interaction of 
Eq. ( |3.152[ |. 
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For weakly interacting pions, the solutions are small perturbations about these values. 
Thus, the rescaled momenta are 



knL 



where Aqn is the small perturbation. 
To first order in Aqn and A 



47r- 



tan^(g„) = nAq^ + OiAq^) 



(3.157) 



(3.158) 



and Eq. (3.1211 yields 



Aq„ 



A 



1 



32^2^ ly^L 



(3.159) 



where Wn'' ~ 2 y + {kn'')'^ is the energy of two free pions with opposite momenta 

of magnitude kn''. Thus, the energy of the corresponding two-pion state, in the presence 
of interactions, is 

H^„ = W^W f I + a'^^— +0{xA . 

" \^ 2 (wrL)3 7 

3.12.4 Two-pion energies from perturbation theory 



(3.160) 



In perturbation theory, the two-pion energy corresponding to pions whose momenta 
would have magnitude kn \ n = 1, ... ,6, in the absence of interactions, can be ex- 
tracted from the tttt — >■ tttt correlation function 



C^.^^^it) = {O„{t)On{0)) c 



where 



E 



{fc„} 



(3.161) 



(3.162) 



is an operator which has overlap with zero-momentum, cubically invariant, two-pio n 
states and t'n is the number of momenta fc„ such that |fc„| = fcj?' (see Eq. (3.150)). 
The sum in Eq. (3.162) is over these momenta, all related by cubic transformations. 
The operator tt{x) has overlap with single pion states. In the limit of large t, the 
contribution of the two-pion states, |7r7rl)y, to the correlation function of Eq. (3.161), 
is 



(3.163) 



where the ellipsis stands for terms which decay more rapidly. The states |7r7r/)v are 
normalized to one. 
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Hs{y) 



O„(0) 



O„(0) 



O„(0) 

(a) (b) 

Fig. 3.5 Diagrams which contribute to C,r7r^7r7r(t) at ©(A") (a) and 0(\) (b). 



A straightforward calculation of the diagrams of Fig. 3.5 using the propagator of 



Eq. (3.154), gives for the correlation function of Eq. (3.161), at ©(A), 



C'ttTT— f TTTT (0 I r 



2 / L3 



1- A- 



i?A((e)^(fc<^'')^ 



L3 



1 



J+o(a2)| 
i^A((fc^"')^(fci°^)^) 



(3.164) 



©(A^ 



i^n \W.nJ " r ■ • 2L3 I (H^^°')3 2 

where we have only retained the contribution which decays exponentially with the rate 



corresponding that of the two-pion state, |7r7rn)y. In Eq. (3.164), the two-pion energy, 



Wm is the same as that obtained from the finite- volume formula of (Liischer, 1991) 



(see Eq. (3.160)) and the regulator contribution is given by 



4(^p + £p) 1 



(3.165) 



with self-explanatory notation. 

Comparison of this result with Eq. (3.163) further gives the matrix element of 
between the vacuum and the two-pion state, |7r7rn)y: 



|(0|O„(0)|7r7rn)y| 



2^ / 



1 + A 



4L3 



1 i?^((e)^(fc(;^')2) 



0(A2) 



(3.166) 
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^ -Hsiy) 



■a 



O„(0) 

Fig. 3.6 Example of a tadpole contribution to C, 



TTTT — ^TVTV 



{t) at 0(A). 



The observant reader will have noticed that we have not taken into account contri- 



butions from diagrams such as the one of Fig. 3.6 which also appear at 0{X). As can be 



verified explicitly, these diagrams amount to a shift of the pion mass by terms which 
are independent of L, up to exponentially small corrections. Since such corrections 
are neglected here, these contributions will affect none of our finite-volume results, 
once the mass has been appropriately renormalized in infinite volume. The details of 
this renormalization are irrelevant here and we assume that the renormalization has 
been adequately performed. Furthermore, the coupling A and the field Tr{x) only get 
renormalized at O(A^), which is beyond the order at which we are working. 



3.12.5 Matching of finite to infinite matrix elements using Eq. (3.146) 



Here we consider a weak transition between the state of a kaon at rest and a two-pion 
state, \mrn)v, in finite volume. The amplitude for this transition is 



M= / v{TTTTn\nw{0,x)\K)v 



The corresponding infinite-volume transition amplitude, T, is given by 
T ^ {7r{p)7T{-p},out\HwmK{6)) . 



(3.167) 



(3.168) 



Again, infinite- volume states are relativistically normalized here. 

To compare with the perturbative results obtained below, we must compute the fac- 



tor rela ting |T| and \M\ in Eq. ( |3.146| to 0(A). Using the expressions in Eqs. ( |3.157[ ) 

0(A2 



(3.159) for qn, and after some algebra, we find 

A 



9n0 (gri) = < 1 



1 

where z„ is the constant given by 



Zn = hm <^ V47rZoo(l;q ) 



— n 



(3.169) 



(3.170) 



Now, using the resuh of Eq. ( |3.156[ ) for 5{k) and Eqs. ( |3.157| )-( [3l59l ) for /c„, we 
find 
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On{tf) 



On{tf) 



Onitf) 





Hs{y) 




K(-U 



K{-ti) K{-ti) 

(a) (b) (c) 

Fig. 3.7 Diagrams which contribute to CK^Trn{ti,tf) at 0{\^) (a) and 0(A) (b,c) 

/ An 



k S'(k )--^^ 



1 — n 



0(V) 



(3.171) 



in Eq. (3.146) 



Combining Eqs. (3.169 1-(3. 171 ), we find for the factor which relates |r| and |M| 



in < q 



1 



dq 
A 1 



2 W}l'^L 





MkW^ 




1.3 






4^ 





(3.172) 



0(A2 



where we have used the expression for fc„ obtained in Eq. (3.157) and (3.159). It should 
be remarked that, unless the pions interact strongly, the size of this factor is essentially 
determined by mismatches in the definitions of T and M and in the normalization of 



states in finite and infinite volume (see Eq. (3.149) and subsequent discussion) 



3.12.6 Matching of finite to infinite matrix elements from perturbation 
theory 



The relevant correlation function here is 

CK^^AU,tf)= [{Onitf)nwio)K{-u,x)) . 

J X 

At 0{X^), it is given by the diagram in Fig. |3.7|a. One trivially obtains, 



2Mk {W^f^Y 



(3.173) 



(3.174) 



At 0(A), the contribution from the diagram in Fig. |3.7[ b gives, using the pion propa- 
gators defined in Eq. (3.154), 



A 



9 e 



2 2Mk 



1 

L3 



Y.^s^^^v'^p)Y X {sAtf - ty- K)Y 
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t 



1 



1 



where 



(3.175) 



(3.176) 



and where we have only kept, in the second hne of Eq. (3.1751, the terms m the 
momentum sum of the firs t line which fall off as e~^" *^ . 

The sum in Eq. (3.176) is restricted to momenta p such that |p| ^ kn'^. To evaluate 
it, we use the asymptotic, large- volume expansion of (Liischer, 1986). Up to terms that 
vanish more rapidly than any power of 1/L, we find 



^(fc(«')^/i(fc(°)) + 



(o)^3 



(3.177) 



where z„ is defined in Eq. (3.170). In Eq. (3.177), I\{kn'') is the infinite- volume con- 
tribution: 



1 



(3.178) 



As discussed following Eq. (3.166), tadpole diagrams such as the one of Fig. 3.7 



solely contribute to the renormalization of the mass of the corresponding leg and do 
not affect our finite- volume expressions. 



Combining Eqs. ( pU74| , ( |3T75| and ( pU77| ), we find, at 0(A), 



K- 



.{U,tf;ki°^)^-g{l--h{k<,^y) 



A 



X <1- 



2^2 VF^")^ 



^i?A((fci'")^(fc^! 



2Mk Wl 

(0)^2^ 



(3.179) 



0(}? 



where, again, W^„, is given by Eq. (3.160). 
3.12.7 Putting it all together 

The contributions of the two-pion states, |7r7rZ)v, to CK^Ti-w(ti,tf) are, in the limit of 
large and tf, 



K- 



rniU,tf) ^e-^^-*'-^'*^ (0|O„(0)|^7r0y v{7nrl\ / nw{0,x)\K)v 



1=0 
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X v{K\K{0)\0) 



(3.180) 



where the ellipsis stands for terms which decay more rapidly. In Eq. (3.180), \K)v is 
a zero- momentum state and, again, all states are normalized to 1. With these normal- 
izations, 



Combining this matrix element with Eqs. (3.179 1-(3. 180) and the result of Eq. (3.166) 
for |(0|C'„(0)|7r7rn)y|, we find 



|M|^i,{l-^/.«»)}i/ 



1 - 



X 1 



4 W/(°'L 



4^ 



(3.182) 



Now, a straightforward evaluation of the infinite-volume analogs of the diagrams 
of Fig. [3Jl yields: 




1 



S \Ep{jP ~ {kil'^y - ie) 



(3.183) 



m 

\M\ 



A 1 



4 WiiU 



47r2 



(3.184) 



which is in perfect agreement with the result of Eq. ( 3.172[ ), predicted by the finite- 
volume formula of Eq. (3.146). In obtaining Eq. (3.184), we have used the fact that 
\T\ = |ReT| + 0(A2). 

3.12.8 K TTTT in finite volume: physical kaon decays 

For illustration, let us suppose that the S'-wave scattering phases Sj, I = 0,2, are 



accurately described by the one-loop formulae of chiral perturbation theory (Gasser 



and Meissner, 1991 Knecht et al., 1995). The two-pion energy spectrum can then be 



calculated in the isospin I channel and in a box of size L where level n — 1 coincides 



with the kaon mass. With this input, the proportionality factor in Eq. (3.146) is easily 



evaluated and one ends up with (cf. Table 3.1) 



|Ao| =44.9x |Mo|, 
1^2 1 =48.7x IM2I, 
IA0M2I = 0.92 X \Mo/M2\ 



(3.185) 
(3.186) 
(3.187) 



As these results show, the large difference between the scattering phases in the two 
isospin channels (about 45° at k — k-^) does not lead to a big variation in the propor- 



tionality factors. In fact, if we set the scattering phases to zero altogether, Eqs. (3.148 )- 
(3.149) give \Ai\ = 47.7 x |M/| for n = 1, which is not far from the results quoted 
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I 


L [fm] 


q 


qd(j)/dq 


kdSi/dk 





5.34 


0.89 


4.70 


1.12 


2 


6.09 


1.02 


6.93 


-0.09 



Table 3.1 Calculation of the proportionality factor in Eq. (3.1461 at the first level crossing 



above. This may be surprising at first sight, since the interactions of the pions in the 
spin and isospin state are quite strong. However, one should take into account the 
fact that the comparison is made for box sizes L which are greater than 5 fm. Hence, 
it is quite plausible that the finite-volume matrix elements already include most of 
the final-state interaction effects. Apart from a purely kinematic factor, only a small 
correction is then required to obtain the infinite-volume matrix elements, from the 
finite- volume ones. 



The proportionality factor in Eq. (3.146) thus appears to be only weakly dependent 
on the final-state interactions. In particular, if the theory is to reproduce the A/ =1/2 
enhancement, the large factor has to come from the ratio of the finite-volume matrix 
elements Mj. In fact, if you carry out this calculation, you should see an 8% enhanced 
AI = 1/2 enhancement! 

To carry out this calculation you will have to address a couple of issues which have 
not yet been discussed here. The first is that, in the absence of twisted boundary con- 
ditions, at least the first excited tttt energy will have to be extracted from the lattice 
calculation. This requires cross-correlator techniques constructed from operators such 



as the one given in Eq. ( 3.162 1 and solving the resulting generalized eigenvalue problem 
(GEVP), as described in (Liischer and Wolff, 1990). The second issue is the one of the 



renormalization of the lattice matrix elements. This is an important, but fairly techni- 
cal problem, which depends sensitively on the fermion discretization used. It is usually 
referred to as the ultraviolet problem, as opposed to the infrared problem which we 
dealt with here, which is associated with the continuation of the theory to Euclidean 
spacetime and the use of a finite volume in numerical simulations. Unfortunately I will 
not have the time to cover the ultraviolet problem here. This problem has been studied 
quite extensively, and I refer you to the original literature, as well as to the lectures of 



Peter (Weisz, 20101 and Tassos (Vladikas, 20101 in this volume. For Wilson fermions. 



which explicitly break chiral symmetry, but retain a full flavor symmetry, the problem 



has been studied in the following series of papers (Bochicchio et ai, 1985 Maiani et al, 



1987 Bernard et al., 1988 Dawson et al., 1998). For twisted-mass QCD, the reference 



is (Frezzotti and Rossi, 2004). When considering domain- wall fermions with a finite 
fifth dimension, one should follow the renormalization set forth for Wilson fermions. 
However, for domain-wall fermions the required subtractions should be significantly 
smaller, since the chiral symmetry breaking should be significantly suppressed com- 
pared to what it is for Wilson fermions. Regarding discretizations which have the 
full, continuum chiral-flavor symmetry at finite lattice spacing (e.g. overlap fermions, 
or domain- wall fermions with a practically infinite fifth dimension), the renormaliza- 
tion is much simplified and will proceed as in the continuum. Finally, for staggered 



fermions, these issues are discussed in ( Sharpe et al., 1987[ Sharpe and Patel, 1994). 
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To conclude, if you wish to be the first particle theorist to unambiguously see the 
A/ = 1/2 rule in K ^ tttt decays and determine e' with controlled errors, I hope that 
you begun working on the problem immediately after the course was given. If not, you 
would do better to hurry because the RBC-UKQCD collaboration is making quick 
progress on these problems (Christ, 2010a Liu, 2010 Sachrajda, 20106). 



4 

Appendix: integral representation for 



Here we derive an integral representation for the zeta function Zoo(l; 9^) of Eq. (3.123 ), 
which is a meromorphic function of , with poles at = n^, n e Z'^. This represen- 
tation is particularly effective for evaluating numerically the kinematic function (l){q) 
that appears in Martin's two-particle momentum quantization formula Eq. (3.121). It 
differs from the one given in Appendix C of ( Liischer, 1991 ). 
The definition of Zoo(s;g'^) is given by Eq. (3.123): 



47r :h„ 

n£2 



n -q ) 



We define 



n -q] 



(4.1) 



where the sum runs over all n £ such that fp > A with A>Keq^. For Res > 0, 



n2>A' 



(4.2) 



with 



A{s;q^) 



Vis) 



n2>A 



rT2<A 



where the second sum runs over all n g such that fp < A. 
To evaluate the sum in Eq. (4.2), we use a Dirac comb: 



^ J d^xf{x)e 



i27Tn-x 



(4.3) 



(4.4) 
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With fin) = e-*" and 



+00 



dx e 



t ^ * 



Eq. (4.4) yields 



Thus 



with 



7]-\3/2 ^2-2 

e *^ 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



Viewed as a function of n, /(s; n) is singular only for n — when Res < 3/2. In 
that case, /'s integrand goes like ^^"^"^/^(l+ig^ +0(^2) when t 0. Thus, Res > 3/2 
yields the half-plane in complex s for which the expression of Eq. (4.71 gives a finite 
result. For such s we can write 



(s; 0) = ^3/2 /• ^s-5/2 l^gt,^ - 1) + ^ 



3/2 
s-^372 ' 



(4.9) 



which is actually well defined for Res > 1/2 and s 7^ 3/2. Thus, for all s in the half 
plane Res > 1/2, we obtain 



r(s) [s-3/2 

-5/2gt,^-4ii I +A(s;q2) 



l^+/^t.-/^(e*^=-l) 



(4.10) 



where the last runs over all n E 1? such that fp 7^ 0. 

Now, for the case s = 1 which is of interest to us here, it is straightforward to 
compute A(s; (7^): 



E 



E 



n? — ^ n? — (P- 



(4.11) 



Using this result in the expression of Eq. (4.101 for Zqq{s; q ) with s = 1, we obtain: 
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^oo(l;g^) = -7r + 



+ 



1^1 



dt 



-{m-q ) 



m — q'^ 



m—1 



dt 



(e*« - 1) 



dt 



(e*' - 1) 



(4.12) 



where counts the ft G such that fp = m (see Eq. (3.150)). Eq. (4.121 is the 
integral representation that we were after. Using this representation, it is straightfor- 
ward to calculate Zqo{1; q^) numerically with good efhciency and high precision, using 
standard integration routines. 
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